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FOREWORD 


This  report  was  prepared  by  Aeronautical  Research  Associates 
of  Princeton,  Inc.,  Princeton,  New  Jersey,  under  Air  Force  Contract 
F33615-69-C-1272.  The  contract  was  initiated  under  Project  No. 
1367*  "Structural  Design  Criteria,"  Task  No.  136702,  "Aerospace 
Vehicle  Structural  Loads  Criteria."  The  work  was  administered 
under  the  direction  of  the  Air  Force  Flight  Dynamics  Laboratory, 
Research  and  Technology  Division,  Air  Force  Systems  Command, 
Wright-Patterson  Air  Force  Base,  Ohio,  Mr.  Paul  L.  Hasty  (FDTE), 
Project  Engineer. 
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Research  Associates  of  Princeton,  Inc.  with  Dr.  John  C.  Houbolt  as 
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This  technical  report  has  been  reviewed  ana  is  approved. 

Gordon  R.  Negaard,  ujor,  USAF 
Chief.  Design  Crite  ia  Branch 
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ABSTRACT 


Recommended  procedures,  based  on  power  spectral  methods,  are 
given  for  the  design  of  aircraft  for  vertical  atmospheric  turbu¬ 
lence  encounter.  Pour  procedures  are  outlined.  The  first  is  a 
direct  preliminary  design  type,  based  on  specified  rigid  body  re¬ 
sults.  The  second  is  more  detailed  and  makes  use  of  composite 
response  values,  computed  either  from  the  rigid  body  results,  or 
by  specific  frequency  response  analysis.  A  third  procedure  is 
based  on  comparison  with  a  previous  sucre- ssful  doci.  In  t  >3 
fourth  procedure  specific  loads.  exceeocS  e  curves  are  derived  in 
accordance  with  assy-rt-jd  missions.  Recommended  gust  intensity 
values,  proportion  of  time  in  turbulence,  and  scale  values  are 
given.  Por  ease  in  application,  detailed  step-by-step  or^cedures 
are  listed  throughout  the  manual. 

The  philosophy  and  basir  for  deriving  the  recommended  de¬ 
sign  values  are  reviewed  in  part.  Special  sections  are  included 
throughout  which  summarize  the  basic  equations  that  apply  to  gust 
analysis,  give  the  means  for  treating  structural  interaction 
effects,  and  show  how  Kxe  problem  of  nonuniform  spanwise  gusts 
is  treated.  Appendices  are  included  which  give  useful  conversion 
and  data  charts,  and  which  review  various  means  for  establishing 
the  frequency  response  function. 
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SECTION  I 


INTRODUCTION 


Under  contract  with  the  Air  Force  Flight  Dynamics  Laboratory, 
Aeronautical  Research  Associates  of  Princeton,  Inc.,  has  been 
developing  procedures  for  designing  aircraft  for  gust  encounter 
based  on  power  spectral  techniques.  Some  of  the  basic  results 
developed  are  presented  in  references  1,  2,  and  3.  A  continuing 
phase  of  effort  has  had  the  objectives  of  further  developing  the 
procedures  and  in  particular  distilling  all  the  work  that  had 
been  accomplished  and  from  this  distillation  to  prepare  a  manual 
for  the  design  of  aircraft  for  vertical  gust  encounter.  This 
report  represents  the  desired  gust  design  manual.  It  should  be 
noted  that  although  effort  was  restricted  to  vertical  gusts  only, 
rcany  of  the  procedures  discussed  herein  can  be  extended  to  the 
cases  of  horizontal  or  combined  vertical  and  horizontal  gust 
encounter.  The  report  thus  sets  the  stage  for  possible  manuals 
on  these  additional  gust  encounter  situations. 

Many  of  the  results  presented  herein  reflect  the  findings  of 
the  past  year  of  effort,  but  notion*  of  many  previous  investiga¬ 
tions  also  enter  prominently,  in  particular,  those  of  references 
2,  4,  and  3.  In  addition,  results  of  many  other  gust  studies 
were  considered  and  incorporated  wherever  possible;  these  studies 
include:  gust  response  studies  related  to  the  present  investiga¬ 
tion,  references  6,  7,  8;  results  of  the  ALLCAT  program,  references 
9-14,  and  many  of  the  reports  dealing  with  the  description  of  the 
atmospheric  turbulence  environment  such  as  references  15-22. 

It  should  be  noted  that,  since  this  report  is  in  the  nature 
of  a  manual,  much  of  the  information  presented  is  given  without 
details  of  development  and  without  elaboration  of  the  logic  by 
which  some  of  the  choices  were  made. 


SECTION  II 


DESCRIPTION  OF  MANUAL 


Problem  treated.-  The  basic  design  problem  that  is  treated 
in  this  manual  is  that  of  an  aircraft  encountering  atmospheric 
turbulence  that  is  considered  to  be  random  in  the  direction  of 
flight  but  which  is  uniform  in  the  spanwise  direction  as  depicted 
in  the  following  sketch 


Purpose  of  manual.-  The  main  purpose  of  the  manual  is  to 
outline  procedures  for  establishing  gust  load  levels  that  are 
appropriate  for  use  in  structural  design  considerations  of  the 
aircraft.  A  secondary  purpose  is  to  give  the  means  for  establish¬ 
ing  the  expected  load  exceedance  curves  for  possible  use  in 
structural  fatigue  analyses  of  the  aircraft  structure. 

Basic  notions  of  analytical  approach . -  The  intent  of  the 
manual  is  to  specify  more  realistic  gust  design  procedures 
through  use  of  the  power  spectral  methods  of  generalized  har¬ 
monic  analysis.  Basic  ingredients  are  essentially  as  follows. 
Modelling  of  the  atmospheric  turbulence  environment  is  made  in 
terms  of  a  spectrum  shape,  a  r.m.s.  gust  intensity  variable 
a ■  ,  a  scale  of  turbulence  L  ,  and  the  parameter  P  depicting 

tne  proportion  of  flight  time  spent  in  turbulence.  The  response 
is  treated  in  the  frequency  plane  instead  of  the  time  plane,  in¬ 
volving  mainly  the  familiar  linear  equation 

♦x(®)  =  |hxH|2<dwH 
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which  relates  the  output  spectra  of  the  response  variable  x  to 
the  input  spectra  of  gust  velocities  w  through  the  frequency- 
response  function  Hx  ,  defined  as  the  response  x  that  results 
from  a  sinusoidal  gust  encounter  of  unit  amplitude.  This  response 
equation  leads  to  thb  two  basic  structural  response  quantities  A 
and  N0  ,  where  A  relates  the  r.m.s.  output  \alue  ax  to  the 
r.m.s.  gust  input  value  through  the  relation 


and  N0  defines  the  number  of  zero  crossings  of  x  with  positive 
slope  per  second.  Load  exceedance  curves  having  the  form 

» -  Mo  f%> 

are  in  turn  developed  for  use  in  establishing  design  loads.  The 
environmental  parameters  P  ,  a ,  and  L  ,  the  structural  re¬ 
sponse  parameters  A  and  N0  ,  and  the  generalized  load  exceed¬ 
ance  curves  are  thus  the  basic  ingredients  of  the  procedures. 

Credibility.-  It  is  felt  appropriate  to  give  the  following 
foreword  on  the  development  of  this  manual  before  giving  a 
synopsis  of  its  contents.  Many  of  the  results  and  procedures 
given  herein  were  established  by  reasonably  sound  mathematical 
analysis.  Because  uncertainty  still  exists  with  respect  to  some 
of  the  quantities  which  characterize  atmospheric  turbulence,  how¬ 
ever,  such  as  the  magnitude  and  distribution  of ' turbulence  in¬ 
tensities  and  scale,  there  are  places  where  specified  values  can¬ 
not  be  regarded  with  certainty.  In  these  places,  recourse  was 
placed  on  reasoning,  judgment,  and  hypothesis. 


Goal.-  An  underlying  goal  in  the  preparation  of  this  manual 
was  to  make  the  procedures  as  simple  as  possible.  Essentially, 
procedures  no  more  Involved  than  the  discrete-gust  design  method 
were  sought.  The  design  approach  evolved  is  basically  as  follows: 
check  the  design  by  the  simplified  procedures  given  (which  are 
analogous  to  the  discrete-gust  design  method);  if  the  design 
passes  the  check,  no  further  work  is  necessary;  if  design  un¬ 
certainties  are  revealed  by  the  simplified  procedures,  then  pro¬ 
ceed  in  the  more  detailed  way  outlined,  but  only  for  those  points 
of  the  structural  design  which  appear  to  be  questionable.  It 
should  be  noted  specifically  that  the  procedures  given  herein  do 
not  depend  on  notions  such  as  probability  of  failure  or  failure 
rate.  There  is  no  need  to  introduce  these  concepts,  and  they  are 

thus  avoided. 
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Contents.-  The  breakdown  and  nature  of  the  contents  of  the 
report  are  essentially  as  follows.  Section  III  lists  the  various 
spectral  relations  that  apply.  Section  IV  defines  the  atmospheric 
input  spectrum  that  i3  recommended.  Section  V  covers  the  practical 
aspects  of  establishing  the  basic  structural  response  quantities 
A  and  N0  .  Section  VI  presents  the  basic  spectral  results  for 
gust  loads  that  apply  to  a  rigid  airplane  with  the  single  degree 
of  freedom  of  vertical  translation;  these  results  form  the  basis 
of  the  simplified  procedures. 

Section  VII  is  the  heart  of  the  report.  It  shows  how,  with 
the  use  of  Section  VI,  to  Identify  the  design  conditions,  how  to 
make  the  preliminary  design  check,  and  then  how  to  proceed  with 
the  cases  that  are  found  to  need  detailed  design  treatment.  Sec¬ 
tion  VIII  outlines  the  procedures  for  establishing  the  basic  load 
exceedance  curves. 

The  intent  of  Section  IX  is  to  show  that  the  procedures  given 
herein  essentially  include  the  considerations  of  discrete-gust 
design.  Section  X  shows  why  limit  load  values  are  chosen  as  the 
base  for  design,  and  why  ultimate  load  considerations  are  avoided. 
The  manner  of  handling  structural  interaction  effects  is  covered 
in  Section  XI.  Although  not  intended  as  part  of  the  presently 
recommended  design  procedure.  Section  XII  is  given  to  show  the 
means  for  handling  the  more  general  gust  encounter  case,  wherein 
the  gusts  are  random  in  the  spanwise  direction  as  well  as  the 
direction  of  flight.  The  section  is  Included  because  It  may  be 
desirable  to  look  Into  this  situation  for  some  of  the  very  large 
aircraft. 

Appendix  A  presents  some  general  data  and  conversion  charts 
that  are  useful  In  spectral  considerations.  Appendix  B  reviews 
the  various  aspects  that  are  involved  in  determining  the  frequency 
response  functions  for  the  aircraft,  the  function  which  is  central 
in  any  spectral  treatment  of  response. 
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SECTION  III 


BASIC  SPECTRAL  RELATIONS 


The  basic  spectral  relations  that  are  of  concern  in  the 
treatment  of  gust  loads  are  presented  in  this  section.  The 
equations  are  presented  in  terms  of  the  frequency  argument  co  , 
but  it  should  be  understood  that  equivalent  expressions  may  also 
be  written  in  terms  of  other  frequency  arguments.  Note  that 
theoretically  a  value  of  infinity  should  be  used  for  the  upper 
limit  in  the  integrals  presented;  for  practical  purpose,  however, 

a  cutoff  frequency  uu  is  used.  The  means  for  establishing  this 

cutoff  frequency  is  given  in  Section  V.  Figure  1,  taken  from 
reference  4,  serves  well  to  illustrate  the  physical  significance 
of  some  of  the  relations  given  in  this  section;  the  notions  of 
this  figure  are  by  now  quite  well  known,  but  they  are  so  basic 
that  reproduction  of  the  figure  here  is  considered  merited. 

Definitions : 

w  represents  the  input  gust  function  w(t) 

x  or  y  represents  the  response  functions  x(t)  and  y(t) 

x  or  y  represents  the  response  functions  x(t)  and  y(t) 

hxH  =  AXH  +  1  bxH 

HyH  =  +  i  By(0)) 

where  Hx  and  H  are  the  frequency  response  functions  for 
the  variables  x  yand  y  due  to  a  sinusoidal  gust  encounter 
of  unit  amplitude;  note  the  A  functions  are  symmetrical 
with  respect  to  the  argument  cd  ,  while  the  B  functions 
are  antisymmetrical.  Various  frequency  arguments  are  as 
follows : 


Frequency  response 


0)  =  2uf 

n  =  “ 

Y 


Hj 


iofix 


(1) 


H •  =  iuH 

y  y 


(2) 
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Basic  input-output  relations 


*x  -  Kl2  K 

(3) 

,  2i„  | 2  . 

0  •  =  a)  H  0 

vx  '  x1 

(*) 

Basic  cross-spectral  relations 

H  H 

0  =0  =  H  H  0  =  — ^  0  =  ■=£  0 

vxy  vyx  x  yvw  Hx  vx  Hy  y 

(5) 

^wx  =  Hx 

(6) 

vxy  x  y 

(7) 

r.m.s.  values  and  correlation  coefficients 

CO 

4  =  /  *wH  d“  (8) 

o 


0) 


0)_ 


Jx=/*x<“>  d“=  /  I 


Hx  *w  du> 


(9) 


ou 


ol 


"  /*y(“)  d“  -  /  i«y|2  K 


da) 


(10) 


0) 


0L 


Paxay 


-  |  f  *xy  (“)  dt"  “  J  Re(RxHy)dw  d“ 


-0). 


0). 


=  /  (A  A  +  B  B  )  0  da) 

J  '  x  y  x  y' 


(11) 
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(12) 


(13) 


(14) 


(15) 


(16) 
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Figure  1.  Input-Output  Relations  for  Gust  Response 


SECTION  IV 


INPUT  SPECTRUM  AND  SHAPE  VARIATIONS 


The  recommended  spectrum  for  the  input  vertical  gust  is  shown 
in  figure  2  and  is  given  by  the  equation 


1  +  |(1.339  LQ)2 

7 - - TJJ76 

1  +  (1.339  IflrJ 


(17) 


Expressions  for  this  spectrum  in  terms  of  other  frequency  argu 
ments  appear  as 


*M(o)  -  $ 
iw(a)  -  l 

^  *y,W 


/  (18) 


The  value  for  the  turbulence  scale  L  has  been  a  very  elusive 
quantity.  Estimates  have  ranged  from  200  ft  to  over  5000  ft  . 
Recent  studies  made  by  the  principal  investigator  of  this  report, 
reference  23 >  aimed  at  deducing  scale  value  by  means  other  than 
those  commonly  used,  have  indicated  that  scale  values  may  be 
smaller  than  those  often  quoted.  On  the  basis  of  these  studies, 
the  value  for  the  effective  scale  of  turbulence  L  is  recommended 
as  L  *  750  ft  for  all  altitudes.  Note  the  scale  is  not  really 
known  to  within  the  accuracy  of  50  ft  ,  or  even  to  within  100  ft  . 
The  choice  is  made  simply  to  make  the  quantity  2L/c  ,  which  is 
frequently  encountered,  expressible  in  a  rounded-off  number  form; 
thus,  the  following  ratio  is  recommended 


2 L 
c 


Useful  information  related  to  the  spectrum  shown  in  figure 
2  is  shown  in  figures  3  through  6.  Figure  3  is  the  type 
of  plot  that  should  be  kept  in  mind  for  all  response  spectra  that 
are  given  herein  and  for  any  that  are  derived  in  design  studies. 
Plots  of  this  type  can  be  quite  significant;  thus,  when  the  spec¬ 
trum  values  are  multiplied  by  the  frequency  and  the  results  are 
plotted  with  a  linear  scale  vertically  and  a  log  scale  horizon¬ 
tally,  then  the  geometric  area  one  sees  visually  under  the  curve 
shows  what  frequency  range  is  contributing  most  to  the  total 
power,  or  mean  square  value.  This  result  is  demonstrated  as  follows 


1.1 


00 

0 

dx 


where  dx  =  du)/o)  or  X  =  log  . 

Figures  4  and  5  give  the  truncated  r.m.s.  and  zero 
crossing  values  of  the  gust  spectrum  as  a  function  of  cutoff  fre¬ 
quency  as  established  by  the  following  equations 

(IA)C 

°c  =  j  *w(W)d(lfl) 


(Lfl), 


(L0)0  = 


/  (Lfi)2  A(LO)d(LO)| 
J0  w 


which  are  similar  to  equations  (8)  and  (16) .  Figure  5  in¬ 
dicates  that  the  "zero  crossings"  value  increases  monotonically 
with  cutoff  frequency.  The  reason  is  that  equation  (17)  is  in¬ 
tended  to  apply  only  in  the  turbulence  inertial  subrange.  A 
leveling  of  the  curve  would  occur  if  a  more  appropriate  equation 
were  used  at  the  high  frequencies  (frequency  components  having 
wavelengths  smaller  than  or.e  or  two  centimeters).  From  a  response 
point  of  view,  however,  the  aerodynamics  and  structural  behavior 
of  the  aircraft  tend  to  wash  out  the  effects  of  high  frequency 
gust  components  even  in  the  inertial  subrange.  Thus,  only  fre¬ 
quencies  below  some  cutoff  frequency  are  of  practical  concern  in 
response  calculations.  Section  V  on  A  and  N  determination 
develops  this  point  specifically. 

Figure  6  gives  a  second  truncated  r.m.s.  value  of  the  gust 
spectrum,  often  used  in  the  experimental  deduction  of  scale  valuq 
as  defined  by 


°t  mJlT 
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Note  that 


2  2  2 
%  =  °c  +  at 


This  section  is  closed  with  some  comment  on  the  variations 
in  shape  of  gust  input  spectra.  In  general,  too  much  attention 
has  been  placed  on  trying  to  fit  various  forms  of  spectral  equa¬ 
tions  to  gust  data.  Often  the  raw  data  are  wrong  cr  unreliable 
to  begin  with,  and  often  the  deduced  spectral  values  are  in  error 
because  of  the  distortions  that  are  produced  by  the  numerical 
analysis  procedure  used.  Usually  there  is  little  or  no  attempt 
to  investigate  how  errors  in  the  data  or  in  the  processing  can 
influence  the  spectral  results,  particularly  at  the  low  frequency 
end.  Thus,  apparent  changes  in  the  spectrum  shape  may,  in 
reality,  not  be  there.  It  should  be  realized  that  considerable 
variations  in  the  apparent  gust  spectrum  shape  can  be  had  even 
using  equation  (17)  as  the  basic  ingredient.  As  the  r.m.s. 
value  of  gust  intensity  varies  from  one  turbulence  encounter  to 
another,  so  too  can  the  scale  be  expected  to  vary.  Thus,  if  a 
distribution  curve  p(L)  is  considered  for  scale  L  ,  see  figure 
7(a),  analogous  to  the  introduction  of  a  probability  density 
distribution  for  ,  then  an  effective  spectrum  shape  may  be 
deduced  using  equation  (17)  as  a  base  according  to  the  equation 


This  deduction  is  recognized  as  being  quite  sijnilar  to  the  deduc¬ 
tion  of  the  familiar  peak  count  relation  e“x'a  from  the  peak 

1  x2 
'  2  2 

count  expression  for  normal  processes  e  °  .  The  following 

example,  taken  from  reference  24,  illustrates  the  use  of  equation 
(19).  Consider  that  p(L)  is  given  by  the  expression 

p (L)  =  .75  6(L  -  400)  +  .25  6 (L  -  2000) 


Effectively,  this  expression  indicates  that  a  normalized  spectra 
having  a  scale  of  400  feet  is  added  to  a  normalized  spectra 


having  a  scale  of  2000  feet,  with  a  weighting  factor  of  .75  on  the 
first  and  .25  on  the  second.  Results  obtained  from  this  expression 
through  use  of  equations  (17)-(19)  are  shown  in  figure  7(b).  The 
effective  spectrum  so  deduced  is  seen  to  contrast  markedly  with  the 
spectrum  given  by  equation  17,  curve  B,  even  though  this  equation 
was  the  basic  building  block  of  the  derived  results.  The  effective 
scale  value  of  the  combined  spectra  is  L  -  800  ft.  Fy  this 
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process  many  different  results  of  the  type  indicated  by  curve 
A  may  be  deduced.  The  results  shown  in  figure  7  thus  indicate 
that  care  must  be  taken  when  atmospheric  spectral  data  are 
interpreted.  Often  it  is  stated  that  the  data  being  presented 
differs  from  other  data;  the  interpretation  is  given,  for  ex¬ 
ample,  that  the  spectra  seem  to  exhibit  a  "mild  knee"  behavior. 
In  reality  the  results  may  simply  be  indicating  that  a  mixture 
of  data  of  several  scales  is  being  presented.  Effectively, 
the  results  are  probably  reflecting  the  fact  that  the  data  are 
not  stationary. 


14 
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Figure  5-  Zero-Crossing  Frequencies  for  Gust  Spectrum 


SECTION  V 


DETERMINATION  OF  A  &  NQ 


Upper  limit  coc.-  Practical  aspects  for  determining  A  and 

by  equations  (15)  and  (16)  are  discussed  in  this  section. 

In  general,  plots  cf  A  and  N0  should  be  made  as  a  function  of 
the  upper  limit  of  integration  o^,.  .  In  some  cases,  no  particular 
difficulty  is  encountered  in  establishing  the  values  of  A  and 
N0  .  Results  may  appear,  for  example,  as  shown  in  figure  8(a); 
the  values  of  A  and  NQ  are  taken  as  the  values  where  the 
functions  have  leveled  off.  In  other  cases,  particularly  de¬ 
pending  on  the  type  of  aerodynamics  used  in  the  establishment  of 
the  frequency  response  functions,  there  does  not  seem  to  be  a 
limiting  value  for  N0  ,  (as  in  the  case  of  figure  5  for  the 
gust  spectral.  The  manner  for  handling  these  cases  is  shown  in 
figure  8(b).  The  top  of  the  knee  of  the  A  curve  is  found; 

N0  is  then  taken  at  the  frequency  value  associated  with  this 
knee.  Essentially,  this  choice  limits  concern  to  only  those 
frequency  components  which  are  significant  in  contributing  to  the 
overall  response.  The  knee  of  the  curve  may  be  tied  in  with 
practical  observations.  Thus,  when  a  record  is  taken  of  some 
response  quantity  such  as  stress,  high  frequency  components  can 
usually  be  seen  in  the  record  trace,  provided  the  instrumentation 
and  recording  system  is  sufficiently  precise.  The  high  frequency 
components  may  be  interpreted  to  be  mainly  a  "hash"  or  noise; 
they  contribute  little  to  the  r.m.s.  value,  but  greatly  affect 
the  N0  count.  The  choice  of  the  top  of  the  knee  has  the  effect 
of  ignoring  or  by-passing  this  high  frequency  "hash"  content  of 
the  record. 

Typical  examples  from  actual  aircraft  applications, 
illustrating  both  of  the  situations  depicted  in  figure  8,  are 
shown  in  figures  9  through  11.  These  practical  applications 
serve  also  to  bring  out  two  facts  of  signlf icance .  It  is  noted 
that  the  value  of  A  is  quite  sensitive  to  the  value  of  scale 
used,  generally  the  higher  the  3cale  the  lower  the  A  value. 

Thus,  A  values  should  not  be  compared  unless  the  same  L 
value  is  used.  In  contrast,  the  value  of  N  is  noted  to  be 
quite  insensitive  to  the  value  of  scale. 
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CUTOFF  FREQUENCY  -  CPS  CUTOFF  FREQUENCY 

Figure  10.  Example  A  and  NQ  Values  Found  for  the  B-58 

Airplane,  5#  Fuel,  50,000  Ft.  (Taken  from  Ref .  6) 
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SECTION  VI 


IC  SINGLE  DEGREE  “OF  “F 


MatrtiMMrtafdliWa 


The  case  of  a  rigid  airplane  with  the  single  degree  of 
freedom  of  vertical  translation  is  treated  in  this  section. 
Response  results  for  eg  vertical  acceleration  in  g-units  are 
established;  these  results  form  the  basis  of  the  preliminary- 
design  consideration.  For  2-dimensional  incompressible  flow* 
the  square  of  the  modulus  of  the  frequency  response  function 
for  this  case  can  be  shown  to  be  (see  reference  25) 


H(k)|2  -  (f^rf  flW  f2<k) 


where 


*!<*)  - 


+  ( + 


f2(k)  =  P 2  +  Q2 


in  which 


2W 

apcgS 


The  function  F  and  G  in  fn  refer  to  the  Theodorsen  func¬ 

tions,  as  used  in  oscillatory  treatment  of  2-dimensional  air¬ 
foils,  while  P  and  Q  in  f2  are  corresponding  type  functions 

for  sinusoidal  gust  encounter.  For  solution  the  function  f2  , 
as  given  in  reference  26,  for  M  -  .2  was  used,  specifically 


p  +1.5*  +  ir  Mk 


2  2 

p  =  1  -  iF 


1  +  4.92k  +  2.o6kc 


M  =  .2 
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The  output  spectra  for  this  case  follows. as 


■Wk>  -  (^ir)2  fi(k)f2(k)*w(*) 


■  W2  <  <*> 

0 

where  ^  =  f-^f^  -vj-  .  From  this  equation  the  r.m.s.  value  for 

aw 

eg  acceleration  in  g-units  is  found  to  be 


°An  ^  K0  °w 


(22) 


where 


Equation  (22)  is  recognized  as  being  analogous  to  the  discrete- 
gust  design  formula.  The  response  parameter  A  for  this  case 
is  seen  to  be 


A 


r 


(23) 


which  also  may  be  expressed  in  the  interesting  form 


a  =xA 

r  eg  n 


(24) 


The  various  results  obtained  for  this  case  are  shown  in  figures 
12  through  18.  Figure  12  shows  the  plots  of  k0h(k)  vs  k  j 
the  curves  of  this  figure  are  significant  because  they  show  the 
ranges  of  k  which  contribute  mainly  to  the  response.  The  k 
values  associated  with  the  maximum  values  of  the  k0n (k)  curves, 
the  location  where  there  is  maximum  contribution  to  the  r.m.s. 
value  of  response,  are  shown  in  figure  13.  The  ky  values 
shown  in  this  figure  are  discussed  subsequently. 
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Figures  14  and  15  show,  respectively,  the  variation  of 

_  —  .  .  ■  J  i  rr 


and  kn  with  cutoff  frequency.  Maximum  values  of  the 

r  _  ^  _  -I.  M  i  .  1  1  _ L  J Zt 1 TT 


.  -  u  -  -  - 

curves  of  figure  14  form  the  gust  alleviation  factor  for 

the  spectral  approach;  these  maximum  values  are  shown  in 
figure  1 6.  For  comparison,  the  gust  factor  K„  as  used  in 
discrete-gust  studies^  reference  27,  is  also  shown.  Figure  17 

gives  the  values  of  —  as  used  in  the  alternate  response 
equation  given  by  equation  (24). 


Figure  15  shows  the  variation  of  kQ  with  cutoff  frequency 
as  established  by  equation  (16).  The  "realistic"  zero-crossing 
values  of  kQ  ,  as  established  from  the  K,*  and  k0  curves  of 
figures  14  and  15,  through  use  of  the  top  of  the  knee  procedure 
described  in  Section  V,  are  shown  in  figure  18.  This  figure 
allows  for  the  quick  evaluation  of  rigid  body  N0  values.  The 
ku  curves  shown  on  the  right  of  figure  13  designate  the  fre¬ 
quency  values  that  were  chosen  for  the  top  of  the  knee.  Es¬ 
sentially  ku  designates  the  upper  value  of  the  range  of  k 
that  is  significant  for  response;  that  is,  there  is  little 
contribution  to  the  r.m.s.  value  for  response  beyond  these 
values.  To  be  systematic  in  selecting  the  top  of  the  knee  for 
this  study,  the  top  of  the  knee  was  specifically  selected  at 
the  point  on  the  K1  curve  where  the  K*  value  was  two  percent 
less  than  the  maximum  value  of  the  curve. 


It  is  of  interest  to  note  that  the  kQ  values  of  figure 
18  are  probably  closely  related  with  the  concept  of  gust  gradient 
distance  that  has  been  used  in  interpreting  gust  encounter  in  a 
discrete  gust  sense;  specifically,  the  average  gust  gradient  that 
has  been  noted  is  felt  to  be  a  reflection  of  the  kQ  values. 

From  kQ  =  gY  ,  and  that  fact  the  ^  is  related  to  wavelength 

by  ^  ,  the  following  relation  of  "equivalent"  gust  gradient 

distance  in  chords  results 


H  _  r 
c  4k 

o 

where  H  refers  to  ^  the  wavelength  X  .  By  this  expression 

the  k0  values  of  figure  18  have  been  converted  to  the 

values  shown  in  figure  19.  It  is  noted  that  for  n  around  20, 
which  is  representative  of  some  of  the  mass  parameters  used  in 

H 

early  gust  studies,  that  —  of  about  10  is  indicated,  in  good 

agreement  with  the  gust  gradient  distances  that  were  noted  in 
the  past  studies.  The  chart  indicates  that  the  average  gust 
gradient  distance  can  be  expected  to  be  greater  for  larger  p. 
values,  and  smaller  for  the  small  pi  values. 
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The  results  presented  in  this  section  are  intended  to  show 
the  character  and  make-up  of  the  acceleration  response  for  the 
rigid  airplane  treated,  and  are  given  for  ready  reference  pur¬ 
poses  in  connection  with  any  response  study  that  may  be  made. 
Equation  (23)  and  figures  16  and  18,  or  alternatively  equation 
(24)  and  figure  17,  are  the  significant  results  of  this  section; 
they  form  the  basis  of  the  preliminary  design  treatment  given  in 
the  subsequent  section. 
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Figure  13.  Frequencies  Associated  with  Peaks  of  the  k0x(k)  Curves 


o 

o 

®o  o 

~  O  o 
ft  tf  CM 


2  to 


~t  d  tt tttzL 


7TTTX _ i_4_  4-jH-i  'I - f  TTi  f  "*  T 

[Hr  r fridti  l;;j± jin  X-ff  {4-1}  ”  * 


4 :i-4-  4—  ^J-i -  lSj  zrtu:xt!:r  if) 

i-.q:3§  ii  VX’-ii 

r  i-H-hrl  IXrr^.j  ...L _ jiP-  -i  -  in  4-+  h  + 


S'q‘14 


mi 

-r  r-|4  4  "ft 

■”1  !"  M  i  T*" 


"" 


H--44- 

M 

P4 


t]  tffitfii  ^ 

*  tifteii  4 


V'\  #  ri4  iV- 
•i  44  -+fir;±t1::M 


-Hit*  h 


iMiipMfli 

•iTX  i  _j±  TfijT  J:j-  4  Iji  •  -ItH:  ;H  ;r 
: j.  hi~  ]  j~- -iti  ‘  *-[£  t  \  \  T  :  •+ 1 1 
t  Tt Jr.  T-  44-’" I"  "TTF”  m-i  Fhx  4-4 
.1  -4  it::]:::  if.;:  |4  U-i  n] . 

i  it  j  4-  -S-4--  i  •  -  i '  f  \  tizm 

flisftHHiPtin' 


-]  -i-i-TTT-  44-  •  • 

1  -:|1  it::']]: 

1  tt  -j  |  it  -J-4-- 

rr tt  iji: -i-t- i 

it  m  P 

J_  -'i  •'  i1  i-  r  ~r' 

■i  •‘i-HT  i’ll'' 


1  M  4 


1,1  •  I  ’ 


1  i,t"  ~T 

■!;  \\\»  '  '  H 

i  i  i  i  !  i  i  I  l  ■  !"■  ■  i  'i  • 


Figure  1^  .  (Cont. 


( C  ont . 


46 


Figure  16.  Gust-Alle viaticn  Factor 
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Figure  18.  Zero-Crossing  Values 


SECTION  VII 


OUTLINE  OF  DESIGN  PROCEDURES 


Nature  of  procedures.-  This  section  outlines  the  gust  design 
procedures  that  are  recommended.  To  keep  the  coverage  as  brief 
as  possible,  very  little  exposition  is  given  in  justifying  the 
numbers  given,  or  in  justifying  the  procedures  -  although  reasons 
should  be  reasonably  clear  for  the  most  part.  To  this  end,  the 
procedures  are  given  mainly  in  step-by-step  fashion.  In  general, 
the  response  is  referred  to  by  the  variable  x  ,  signifying  any 
response  quantity  of  concern,  such  as  acceleration,  a  bending 
moment,  or  a  stress.  It  should  be  kept  in  mind  that  the  response 
values  established  are  incremental  values  due  to  gust  encounter, 
and  that  these  must  be  superimposed  on  the  1-g  level  flight 
values  to  obtain  the  total  response.  Note  also  that  in  this 
section  the  response  refers  to  an  individual  or  specific  load 
variable  at  a  particular  point  on  the  aircraft  structure.  The 
means  for  handling  the  situation  of  combined  loads  or  stresses 
is  shown  later  in  Section  XI. 

A  flow  chart  showing  the  basic  sequence  that  is  involved  in 
applying  the  design  procedures  is  shown  in  figure  20.  Study  of 
this  chart  may  help  in  understanding  the  procedural  aspects  that 
follow.  Essentially,  the  procedure  is  as  follows.  It  is  assumed 
that  the  aircraft  components  have  been  sized  or  that  the  load 
values  have  been  established  in  accordance  with  some  prescribed 
maneuver  load  factor.  Loads  or  stresses  due  to  a  1-g  level  flight 
condition  are  in  turn  established  for  each  of  several  flight 
conditions  that  are  likely.  A  preliminary  check  for  gust  en¬ 
counter  is  then  made  by  fairly  simple  procedures,  based  on  the 
rigid-body  results  of  the  previous  section  and  the  1-g  flight 
loads.  If  the  design  fails  this  preliminary  check,  evaluation 
proceeds  into  a  phase  of  intermediate  detail.  This  second  phase 
still  makes  use  of  the  rigid-body  results  of  the  previous  section, 
but  considers  the  various  possible  flight  conditions  in  combina¬ 
tion.  Should  the  design  fail  this  phase,  design  proceeds  to  a 
more  detailed  study  phase  where  explicit  evaluation  is  made  of 
some  of  the  frequency  response  functions.  This  is  the  only  phase 
where  detailed  evaluation  of  frequency  response  functions  is  re¬ 
quired;  the  necessity  of  having  to  enter  this  phase  will  usually 
be  remote.  Note  the  use  of  the  level  flight  load  values  x-j„  in 
the  preliminary  and  intermediate  phases  is  a  key  notion;  "  6 
essentially,  the  sizing  of  structural  members  in  a  direct  fashion, 
rather  than  by  an  iterative  technique,  is  implied. 

The  design  borders  and  conditions  presented  in  this  section 
are  based  to  a  large  extent  on 

a)  discrete-gust  design  values 

b)  j  and  NQ  values  as  obtained  from  airplane  computation 

studies,  references  6-8,  and  the  related  analysis  given 
in  reference  2 
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c)  limit  load  exceedance  values  given  in  reference  1 6 

d)  results  given  in  reference  28 

e)  implications  obtained  by  carefully  comparing  and  inter¬ 
preting  the  results  of  past  gust  studies. 


To  derive  procedures  in  specific  detail,  the  following  equations 
for  load  exceedance  were  adopted  as  a  base  for  establishing  the 
main  results  of  this  section 


If  -  PN  e 
o 


(25) 


n  =  PTN  e 
o 


(26) 


These  equations  are  assumed  to  apply  in  the  limit  load  region  of 
response.  From  equation  (25)  the  following  interesting  alternate 
forms  may  be  derived 


x_ 

a 


PVk 


Nh  =  3600Pttoe  X  -  3222  -^S.  e 


x_ 

a. 


N  = 
m 


N„  = 


PN  ' 
5280  -y  e 

X 

—  Pk  e  0x 
7T  O 


x  =  5280  o  -x 

7T  C 


(27) 

(28) 

(29) 


where  Nh  denotes  the  number  of  upward  crossings  per  hour  of 
travel,  ^m  denotes  the  number  of  upward  crossings  per  mile  of 
travel,  and  Nc  denotes  the  number  of  upward  crossings  per  chord 
of  travel.  For  various  reasons,  not  discussed  here,  representation 
of  load  exceedance  based  on  Nh  is  preferred.  It  is  interesting 
to  note,  however,  that  of  the  several  forms,  the  equation  for  Nc 
seems  to  be  the  most  "natural"  form,  in  the  sense  that  it  is  the 
only  form  which  is  completely  nondimens ional,  and  because  it  in¬ 
volves  the  least  number  of  parameters. 

In  addition  to  equations  (25)  and  (26),  the  basic  gust  de¬ 
sign  condition  assumed  in  this  development  is  that  the  number  of 
limit  load  exceedances  shall  be  10  or  less  in  30,000  hours  of 
flight.  The  choice  of  10  limit  load  exceedances  implies  5  for 
positive  loads,  5  for  negative  loads.  Based  on  these  values,  the 
limit  load  exceedance  rate  in  number  per  second  for  satisfactory 
gust  design  is  given  by 


52 


(30) 


N  < 


^  30,000  X  3600 


=  4.62  X 


With  this  expression,  and  the  fact  that  ov  =  AaTI 
may  be  written  X  W 


4.62  x  10"8  > 


PN  e 
o 


equation  (25) 


(31) 


This  equation  forms  the  basis  of  the  design  procedures. 

For  the  preliminary  design  phase,  equation  (31)  is  applied 
to  specific  flight  conditions  and  altitudes.  Design  is  judged 
satisfactory  if  this  equation  is  satisfied  by  the  values  of 

and  N0  at  each  of  these  chosen  flight  conditions.  This  applica¬ 
tion  implies  the  specification  of  the  values  of  P  and  a™  . 
Recommended  values  of  P  and  crw  considered  appropriate  for  de¬ 
sign,  which  were  arrived  at  after  much  deliberation  and  reasoning, 
are  shown  in  figures  21  and  22.  With  these  P  and  ctw  values, 

equation  (31)  yields  the  N0  vs.  j  preliminary  design  borders 

shown  in  figure  23.  The  use  of  this  figure  implies  the  following 
concept.  Consider  that  the  airplane  is  to  be  flown  wholly  at  a 
given  altitude  for  the  entire  lifetime  expected  of  the  aircraft. 

If  the  design  appears  satisfactory  for  all  such  altitudes,  then 
certainly  the  design  should  be  adequate  regardless  of  what  combi¬ 
nations  of  altitudes  are  experienced.  Detailed  use  of  figure  23 
will  be  outlined  subsequently. 

For  the  intermediate  and  detailed  design  phases,  equation 
(31)  is  also  used,  but  care  is  taken  to  combine  the  results  from 
various  flight  conditions  in  a  systematic  manner.  Details  of  the 
processes  will  also  be  outlined  subsequently. 

It  is  of  interest  to  note  that  reference  16  indicates  the 
number  of  limit  load  exceedance  per  mile  of  flight  of  %  =  10“°  . 
This  value,  and  the  chosen  values  of  n  =  10  ,  T  =  30,000  hrs. 
Implies  an  average  flight  speed  given  by  the  relation 

_  _  VT  „ 
n  ~  Nm 


10  =  2Q^QQQ^36o°  10” 6  v 

This  relation  yields  V  =  488  fps  ,  which  appears  representative 
of  operational  practice.  The  values  of  n  =  10  ,  T  =  30,000  hrs, 
and  Nm  =  lo”°  thus  appear  consistent. 
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Other  values  considered  representative  of  a  sound  gust 
design  condition,  and  which  are  implied  by  the  chosen  design 

numbers  are  as  follows.  Suppose  ^  is  taken  as  5^5  at  sea 

level  where  aw  =  6  (figure  2i) ,  this  yields  x  -  x 


a. 


*  9 


From  figure  23,  for  h  =  0  and 
This  N, 

as  follows: 


S=  54 


N 


x  Aaw 
is  found  to  be 


Nn  =  .75  ops. 

O)  c 

yield  k„  = 


o  value,  and  the  value  of 


kQ  =  .0966 


‘o  -  2V 

k  =  .0724  for  c  =  15  .  These  k 
o  o 

representative  of  practical  situations, 
the  various  numbers  chosen  or  derived  in  this  section  leads 
to  the  following  combination  which  is  considered  representa¬ 
tive  of  a  sound  gust  design  situation 


V  =  488  , 
for  c  =  20  , 
values  are  also  quite 
Thus,  a  summary  of 


n 

T 


=  5 


30,000  hrs 

io'6 

488  fps 
.0005 


SL  = 

m 

V  = 

P  = 

x 

N  =  .75  cps 


These  numbers  are  self-consistent  and  satisfy  equations  (25) 
through  (28). 

Preliminary  design  considerations .-  A  preliminary  -  and 
perhaps  final  -  design  check  proceeds'  as  follows. 

1,  List  various  possible  conditions  of  flight.  Involving 
altitude,  speed,  weight  and  weight  distribution,  as 
depicted  in  the  following  sketch 
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2.  Select  points  throughout  the  structure  that  are 
suspected  of  being  critical  locations. 

3.  Establish  values  of  the  level  flight  load  or  stresses, 

x1_  ,  at  these  locations  for  the  various  flight 

conditions. 

4.  Establish  the  values  of  Ar  for  each  of  the  chosen 
flight  conditions  by  equation  (23)  or  (24),  using 

the  values  of  K,  or  ,  as  given  in  figures  1 6 

li, 

and  17.  With  these  A  's  establish  An  values  by 

r 

the  expression 

An  =  Ar9aw  (32) 

where  a  is  taken  from  figure  21.  These  An  values 
are  verywsimilar  to  the  values  that  would  be  established 
by  the  discrete-gust  design  approach.  The  main  purpose 
of  this  step  is  to  isolate  those  flight  conditions  which 
are  most  likely  to  produce  the  largest  gust  loads.  Thus 
if  An  is  found  to  exceed  the  value  of  the  design 
maneuver  acceleration  value  minus  one,  such  as  the 
Arn.  =1.5  factor  often  used  in  the  past,  then  the  flight 
condition  yielding  this  An  will  probably  also  be  found 
critical  in  the  subsequent  steps.  If  all  An's  are  less 
than  the  maneuver  acceleration  minus  one  value,  gust  de¬ 
sign  is  probably  not  a  problem. 

5.  With  the  Xy_g  values  of  step  3  and  the  Ar's  of  step  4 
determine  the  values  of  A  by  the  relation 

A=Arxl-g  (33) 

For  aircraft  with  large  flexible  swept  wings,  the  slope 
of  the  lift  curve  used  in  step  4  should  be  that  for  the 
flexible  airplane,  both  in  determining  p,  and  Ar,  so 
‘that  the  load  effects  due  to  wing  bending  are  approxi¬ 
mately  taken  into  account.  Note  the  intent  of  consider¬ 
ing  various  flight  conditions  is  to  find  the  condition 
which  leads  to  the  largest  value  of  A  at  each  altitude 
for  the  particular  response  quantity  of  concern.  Note 
also,  that  there  is  no  explicit  determination  of  the 
frequency  response  function  or  output  spectra  in  this 
preliminary  check  sequence;  results  for  the  rigid-body 
case  of  Section  VI  are  used  directly,  much  as  in  the 
discrete-gust  design  technique. 

6.  For  each  flight  condition  estimate  the  value  of  N 

using  figure  18.  0 


55 


7. 


For  each  structural  point  being  checked  and  each 
altitude *  take  the  largest  values  of  A  found, 
multiply  by  1.1  and  divide  the  resulting  value  into 

XL  “  xl-g  j  where  x^  is  the  limit  load  value  for 
the  response  quantity;  specifically,  the  ratio 


XL  -  *1-1 


is  found. 


becomes  the  effective 


With  equation  (33)  this  ratio 
t  value  as  follows 


Note,  the  factor  1.1  is  introduced  as  a  means  for 
approximately  taking  into  account  the  amplification 
effects  due  to  flexibility.  The  factor  is  a  rough 
average  value;  if  judgment  or  some  previous  results 
indicates  that  a  different  factor  may  be  in  order, 
the  number  may  be  adjusted  upward  or  downward 
accordingly. 

8.  Enter  figure  23  with  the  NQ  and  £  values  es¬ 
tablished  in  steps  6  and  7  and  compare  each  point 
with  the  appropriate  altitude  curve.  Decide  action 
according  to  the  following  three  rules: 

Rule  (a):  If  the  point  lies  well  to  the  right 
of  the  curve,  the  design  may  be  con¬ 
sidered  satisfactory  from  a  strength 
point  of  view,  and  no  further  cheek 
need  be  made. 

Rule  (b):  If  the  point  lies  to  the  left  of  the 

curve,  the  design  is  probably  unsatis¬ 
factory,  and  needs  more  detailed  study, 
as  outlined  in  detailed  design  phase. 

Rule  (c):  If  the  point  lies  close  to  the  curve, 
on  either  side,  the  design  may  or  may 
not  be  satisfactory,  and  also  needs 
more  detailed  study. 

9.  If  either  rule  (b)  or  (c)  of  step  8  is  encountered, 
proceed  to  the  detailed  design  consideration. 

It  should  be  noted  that  the  underlying  thought  behind 
rule  (a)  is  that  if  the  airplane  is  considered  to  be  flown 
continually  in  the  condition  that  leads  to  the  most  severe 
loads  at  a  given  point;  and  the  structure  is  shown  to  be  ade¬ 
quate  for  this  situation,  then  certainly  the  structure  should 
be  satisfactory  for  all  other  flight  conditions. 


Detailed  design  considerations.-  Procedures  for  examining 
the  design  in  a  more  detailed  sense  than  offered  by  the  pre¬ 
liminary  design  considerations  are  listed  here.  The  procedures 
have  some  aspects  of  a  mission  consideration  but  are  somewhat 
easier  because  explicit  evaluation  of  load  exceedance  curves  is 
avoided.  Three  procedures  are  given*  each  of  different  complexity 
and  detail.  The  first  is  a  composite  approach  based  upon  c-g 
acceleration  values.  The  second  is  in  terms  of  limit  load  ex¬ 
ceedance  rates  but  is  based  on  the  rigid-body  results.  These  two 
are  of  intermediate  detail.  The  third  is  also  in  terms  of  limit 
load  exceedance  rates  but  involves  detailed  evaluation  of  the 
frequency  response  functions.  All  three  make  use  of  equations 
(25)  and  (30). 


Composite  approach  based  on  c-g  acceleration: 


This  approach  may  be  used  if  rule  (b)  or  (c)  of  step  8  of 
the  preliminary  design  procedure  applies.  It  is  not  a  basic 
design  method,  but  rather  serves  as  a  check  as  to  whether  a 
more  detailed  design  treatment  appears  necessary.  It  is  given 
because  it  is  quite  simple  to  apply,  and  because  it  does  not 

require  the  use  of  the  x.  load  values. 

j-”S 

1.  Lay  out  the  most  probable  mission  (or  missions)  of  the 
airplane,  as  depicted  in  figure  24.  Usually  missions 
with  a  cruise  altitude  of  about  30,000  ft  will  lead  to 
the  most  severe  gust  load  conditions  (unless  of  course 
the  mission  is  a  low  altitude  one  at  very  high  speeds). 


2.  Establish  the  values  of  Ar  and  N0  for  each  segment 
of  the  mission;  use  either  equation  (23)  or  (24)  and 
figure  1 6  or  17  and  l8j  also  establish  the  value  of 
a  for  each  segment  by  the  relation 

A 


3. 


where  a  is  taken  from  figure  21. 
w 

Determine  composite  values  of  the  parameters  a 
PNq  through  means  of  the  following  equations:  x 


(35) 

and 


T  =  ^  +  T2  +  T3  +  ... 


PT 
2 

c  = 


P1T1  +  P2T2  +  P3T3  +  • • • 


x 


...) 


PN  = 


■  W  (Wxj  +  p2Vx2  +  p3T3°x3 

=  J  (P- Tn N  +  P0T0N_  +  P0T0N  +  ...) 


(36) 

(37) 

(38) 

(39) 


where  the  subscripts  denote  the  segments  and  where  the 
P  values  are  given  by  figure  22.  These  equations 
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represent  simply  the  overall  parameters  that  are 
obtained  when  all  the  segments  are  considered  in 
combination. 

4.  From  equation  (38)  establish  the  effective  value  of 
—  by  the  relation 

c 

X  Vv. 

x_  _ 

°x'^X 

where  An^  is  the  design  maneuver  load  factor  minus 
one.  Enter  figure  25  with  this  ratio  and  the  value  of 
PNq  of  equation  (39) » 

5.  If  the  point  lies  well  to  the  right  of  the  curve,  the 
design  may  be  considered  satisfactory;  if  the  point 
lies  near  or  to  the  left  of  the  curve,  more  detailed 
examination  is  needed  as  in  the  following  two  pro¬ 
cedures.  The  point  shown  in  figure  25  applies  to  an 
example  case  that  is  treated  later. 


Approach  based  on  rate  of  limit  load  exceedances; 


This  approach  is  also  used  if  rule  (b)  or  (c)  of  step  8 
in  the  preliminary  design  section  applies,  or  if  the  check 
given  by  the  preceding  composite  approach  failed.  The  approach 
is  based  on  the  rigid-body  results  of  Section  VI,  and  uses  the 
xl_g  values. 

1.  Proceed  as  in  steps  1  and  2  of  the  preceding  composite 
approach,  but  also  establish  the  x-  values.  From 
the  A  and  Xj~  Values  obtain  tnefe  ox  value  for 
each  segment  by  tne  relation 


cr  =  Ax,  _  cr 
x  r  1-g  w 


(4o) 


where  a  is  taken  from  figure  21  at  the  appropriate 
altitude^ 


2.  Determine  the  rate  of  limit  load  exceedances  by  the 
equation 


N  =  P 


1  T" 


N 


- 

e  +  P 


2  T 


N 


-  a2 
e  + 


where 


OS  a 
n  n 


fc£r-‘] 

6  n 


(41) 
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and  where  the  subscripts  correspond  to  the  respective 
segments. 


3.  If  the  value  of  N  as  obtained  in  step  2  is  signifi¬ 
cantly  lower  than  the  4.62  x  10"°  value  given  by 
equation  (30),  the  design  may  be  considered  satis¬ 
factory.  If  the  value  of  N  is  near  this  value  more 
detailed  consideration  is  needed  as  in  the  detailed 
approach  given  next.  If  N  is  above  the  value  of 
4.62  x  10"®  3  the  design  is  probably  unsatisfactory 
and  should  be  also  examined  in  more  detail. 

Approach  based  on  rate  of  limit  load  exceedances 
(Detailed  evaluation): 

If  the  design  has  failed  all  previous  three  approaches 
(the  preliminary  check,  the  composite  approach,  and  the  rate 
of  limit  load  exceedances  approach),  then  detailed  design 
consideration  should  be  given  according  to  the  following: 

1.  Proceed  as  in  the  previous  approach  dealing  with 
rate  of  load  exceedances,  except  establish  the 
values  of  A  and  N0  in  accordance  with  equations 
(15)  and  (16).  The  previous  approach  using  the  rigid- 
body  results  will  identify  those  segments  which  con¬ 
tribute  most  significantly  to  the  N  value.  Usually, 
one  or  two  segments  will  contribute  N  values  which 
far  overshadow  the  values  of  N  that  are  due  to  the 
other  segments.  Detailed  evaluation  of  A  and  N 
need  only  be  made  for  these  more  critical  segments? 
results  of  the  preceding  section  may  be  used  for  the 
segments  which  contribute  in  a  minor  way  tc  N  , 

Note:  this  is  the  first  place  in  the  procedures  out¬ 
lined  so  far  where  it  is  necessary  to  evaluate  ex¬ 
plicitly  some  of  the  frequency  response  functions, 
and  in  turn  the  associated  values  of  A  and  N  ; 

the  Intent,  however,  is  to  keep  the  number  of  sfich 
evaluations  to  a  minimum.  For  the  A  values  that 
are  evaluated  explicitly,  it  is  not  necessary  to  use 
the  x-  values,  since  the  correct  units  will  in- 
herentxysbe  included. 

2.  Evaluate  N  by  means  of  equation  (4l).  Note:  for 

the  segments  for  which  A  and  N  have  been 
evaluated  in  detail,  use  0 

_  *L  -  Xl-g 
n  A 

n  wn 

That  is,  the  factor  1.1  is  not  included  due  to  the 
fact  that  amplification  effects  due  to  flexibility 
will  be  inherently  present  in  the  A  values. 


59 


3.  If  the  value  of  N  as  obtained  in  step  2  is  equal 
to  or  less  than  4.62  x  10“®  ,  the  design  may  be 
judged  adequate.  If  N  is  greater  than  this  value, 
strengthening  of  the  structure  would  be  indicated. 


Design  using  a  comparison  approach.-  Another  design  approach 
that  may  be  used  involves  the  concept  of  comparing  the  design 
under  consideration  with  a  previous  design  which  has  ’'proven”  it¬ 
self  through  years  of  successful  operation.  This  comparison 
approach  would  seem  to  be  particularly  attractive  if  the  new  de¬ 
sign  is  similar  in  configuration  to  a  previous  model,  or  repre¬ 
sents  a  stretched  version.  The  comparison  approach  is  ac¬ 
complished  through  use  of  figure  2b,  which  is  derivable  directly 
from  equation  (26).  The  parameters  that  are- used  to  enter  this 
figure  are  evaluated  in  accordance  with  the  composite  values, 
equations  (36)  -  (39)-  If  the  point  established  by  entering  the 
figure  with  the  computed  ratio  values  falls  in  the  "safe”  region, 
then  the  new  design  may  be  considered  satisfactory;  if  the  point 
falls  on  the  "unsafe"  side,  then  more  detailed  study  of  the  new 
design  is  indicated.  Note  that  even  though  the  concept  of  com¬ 
paring  one  aircraft  to  another  is  involved,  this  approach  also 
serves  well  even  for  the  same  airplane,  using  the  idea  of  com¬ 
paring  in  a  relative  sense  one  point  of  the  structure  with 
another  point  of  the  same  structure.  In  this  case,  the  PT  and 
c  values  would  simply  cancel.  The  approach  js  thus  useful  in 
establishing  which  point  of  the  structure  seems  to  be  the  most 
critical  in  relation  to  other  points. 

Examples.-  Some  examples  are  given  to  close  out  this 
section  to"  "illustrate  the  gust  design  procedures. 


Example  1: 


Consider  an  airplane  with  characteristics  as  follows: 


b  =  140  ft 

c  =  20  ft 
m 

S  =  2800  ft2 
a  =  4.88 

W£  =  128,000  lbs  (Empty  wt) 

Wp  =  32,000  (Payload  wt) 

Wp  =  140,000  (Fuel  wt) 

W&  =  300,000  (Gross  wt) 

AnL  =1.5  (Maneuver  load  factor  =  2.5) 


Design  considerations  for  wing-root  bending  moment  will  be 
illustrated. 


For  the  preliminary  design  check,  three  flight  conditions 
are  considered,  as  listed  in  the  following  table: 


h,ft 

W,lbs 

V,fps  (B.M.),, ft-lbs 

- 

(1) 

0 

300,000 

422  2.02  x  10° 

r 

(2) 

30, 000 

300,000 

800  2.02  x  io° 

r 

(3) 

30,000 

180,000 

,  800  1.90  x  10 

where  B.M.  denotes  wing-root  bending  moment.  Condition  (1) 
also  applies  to  maneuver  load  design  and,  therefore,  with  a 

maneuver  load  factor  of  2.5,  the  maneuver  design  bending 
moment  is 

xL  =  2.5  x  2.02  x  106  =  5.05  x  106  ft-lbs 

K 

The  various  values  of  ^  ,  --&■  ,  Ar  that  apply  to  each  of 

these  flight  conditions,  together  with  the  values  of  An  as 
given  by  equation  (32),  the  j  values  as  given  by  equation 

(34),  and  the  N0  values  as  determined  from  figure  18,  are 
listed  in  the  following  table: 


M- 

M- 

Ar 

aw 

An 

X 

A 

N^ 

0 

(1) 

28.6 

.0225 

.0147 

6.0 

.794 

92.7 

.517 

(2) 

76.6 

.0103 

.0125 

9.6 

1.08 

109.2 

.725 

(3) 

46.0 

.0158 

.0196 

9.6 

1.69 

76.8 

.828 

where  a  value  of  ~  =  =  75  was  used.  The  first  point  to 

notice  is  the  An  =  I.69  load  factor  for  flight  condition  (3); 
this  value  is  in  excess  of  the  AnL  =  1.5  value.  Flight  condi¬ 
tions  which  approach  condition  (3)*  which  corresponds  to  a 
lightly  loaded  airplane  (l4  percent  fuel)  at  high  altitude  at 
high  speed,  may  thus  be  suspected  as  possible  critical  gust 
loads  producers. 

A  plot  of  the  ^  and  NQ  values  and  the  appropriate 

design  borders  of  figure  23  are  shown  in  figure  27.  Conditions 
(1)  and  (2)  are  quite  "safe"  but  condition  (3)  is  noted  to 
appear  as  "unsafe."  More  detailed  examination  is  thus  indicated 
this  examination  is  indicated  in  the  following  section. 


Detailed  design  check: 


The  mission  chosen  for  study  is  shown  in  figure  28;  a 
take-off  gross  weight  of  300,000  lbs,  involving  a  full  payload 
and  full  fuel  condition,  is  considered.  Near  the  top  of  figure 
28  are  listed  the  operational  parameters,  and  the  1-g  root¬ 
bending  moment.  Next,  values  of  a  and  P  as  obtained  from 
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figures  21  and  22  are  listed.  The  pertinent  response  para¬ 
meters  (p,  jx  ...  N0)  that  are  evaluated  by  the  various 
equations  and  charts  of  this  report  are  then  given. 

The  composite  values  of  ax  and  PN0  as  obtained  from 
equations  (38)  and  (39)  (using  ax  as  obtained  from  equation 
(35))  are  as  follows: 

°x  =  tF  [l5(.087)2  +  15(.106)2  +  120(  .1324) 2  +  120(  .1535) 2 

+  15(.l87)2  +  15(.ll8)2] 

=  (.11  +  .17  +  2.1  +  2.84  +  .53  +  .27) 

=  .0198 

°x  =  *141 

pNo  =  (15  x  .48  +  15  x  .60  +  120  x  .75  +  120  x  .78 

+  15  x  .87  +  15  X  .55) 

-  .000368 

The  values  of  .15  and  300  appearing  in  these  equations  are  the 

PT  and  T  values  as  evaluated  in  figure  28.  The  value  of  — - 
is  thus  found  to  be  x 


9.7 


This  value  and  the  value  of  PN0  plot  as  shown  in  figure  25$ 
since  the  point  is  to  the  right  of  the  design  border,  the  de¬ 
sign  would  appear  to  be  safe. 


A  more  definitive  check  is  afforded  by  the  approach  in 
terms  of  rate  of  limit  load  exceedances.  For  this  example 
equation  (41)  is  found  to  be 

N  «  3§§^(l5  x  .48  e"15,7  +  15  x  .60  e”12'86  +  120  x  .75  e"10,57 

+  120  x  .78  e“9,42  +  15  x  .87  e’7*86  +  15  x  .55  e"12,54) 

=  (*01  +  *19  +  18’7  +  6l'5  +  41  +  -24> 


-  2.51  x  10 


-8 
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The  o^j  values  (the  exponential  values  of  e)  in  this 
equation  are  found  from  the  ax  values  given  by  equation  40, 
the  (B.M. )^_g  values,  and  the  xL  value  of  5. 05  x  10°  ft-lbs. 
The  value  of&8l00  appearing  in  the  evaluation  is  the  result  of 
removing  a  factor  e"°  from  all  the  exponential  terms  in 
parentheses,  so  as  to  simplify  the  evaluation  of  each  term. 

The  value  for  N  that  is  found  is  noted  to  be  well  below  the 
design  value  of  4.62  x  10"°  given  by  equation  (30).  The  de¬ 
sign  is  therefore  judged  safe.  It  is  to  be  noted  that  the 
equation  for  N  shows  that  the  4th  and  5th  segments  are  the 
largest  contributors  to  the  load  exceedance  rateQ  If  N  were 
found  to  be  larger  than  the  value  of  4.62  x  1CT°  ,  then  evalua¬ 
tion  would  proceed  by  evaluating  the  frequency  response  function 
for  bending  moment,  and  the  A  and  N0  values,  for  only  the 
4th  and  5th  segments.  This  observation  thus  brings  out  the  fact 
that  the  design  procedures  given  herein  are  intended  to  keep  the 
amount  of  detailed  work  to  a  minimum. 

Example  2: 

In  this  example  let  all  the  conditions  be  the  same  as  in 
example  1,  but  with  changes  in  segments  4  and  5  as  follows. 
Assume  that  previous  work  with  the  airplane  design  indicated 
that  the  1.1  factor  chosen  for  flexible  body  amplification 
effects  was  too  small.  Thus  specific  evaluation  of  the  fre¬ 
quency  response  functions  and  the  A  and  N0  values  appeared 
necessary.  Suppose  that  these  evaluations  were  made  and  that 
values  of  04  =  8.63  and  etc  =  7.2  were  found  as  a  result 
(contrasted  to  04  =  9.42  ana  cL  =  7.86  of  example  1).  These 
more  refined  a  values  indicated  in  effect,  that  the  amplifi¬ 
cation  effects  due  to  flexibility  were  indeed  greater,  as 
suspected;  specifically  a  factor  of  around  1.2  instead  of  1.1 
is  indicated.  With  these  04  and  otc  values,  and  all  other 
values  as  in  example  1,  the  equation  for  N  becomes 


N  =  366’x°Bli6o  (-01  +  -19  +  l8-7  +  *35  +  79  +  .24) 

=  4.8  x  10’8 

This  value  is  greater  than  the  design  value  of  4.62  x  l(f8 
given  by  equation  (30),  The  value  is  only  4#  greater  and 
thus  the  design  appears  to  be  only  slightly  inadequate.  If 
strict  adherence  to  the  design  number  is  maintained,  however, 
some  strengthening  cf  the  structure  due  to  gust  encounter 
effects  is  indicated. 
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If  fall  marginally 


Figure  20.  Sequence  of  Gust  Design  Procedure 


Figure  22.  Recommended  P  Values  for  Composite  Approaches 
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Figure  25.  Design  Chart  Using  Composite  Value 


Figure  27.  Preliminary  Design  Results  for  Example  Problem 


Segment 

J. 

3 

4 

3 

0 

T 

15 

15 

120 

ro 

0 

15 

15 

W 

297.5 

292.5 

270 

230 

207.5 

202.5 

V 

420 

600 

800 

800 

800 

420 

(B.M. )1_g 

2.02 

2.02 

1.99 

1.95 

1.93 

1.92 

ow(fig.  21) 

6.9 

8.7 

9.6 

9.6 

8.7 

6.9 

P(fig.  22) 

.0005 

.0005 

.0005 

.0005 

.0005 

.0005 

P 

.0019 

.00116 

.00889 

.00889 

.00116 

.0019 

H 

35.5 

57.2 

69 

58.7 

4o,6 

24.2 

h 

.0193 

.0131 

.0111 

.0129 

.0173 

.0262 

M- 

Ar 

.0126 

.0122 

.0138 

.0160 

.0215 

.0171 

ax(eq.  35) 

.087 

.106 

.1324 

.1535 

.187 

.118 

ax(eq.  40) 

.1756 

.214 

.264 

.299 

.361 

.227 

k 

o 

.072 

.063 

.059 

.061 

.068 

.082 

No 

.48 

.60 

.75 

.78 

.87 

.55 

T  =  15  +  15  +  120  +  120  +  15  +  15  •=  300 
PT  =  .0005  x  300  =  .15 


min 

x  1000  lbs 
fps 

x  106  ft-lbs 


106  ft -lie 


Figure  28.  Mission  Characteristics  of  Example  Problem 


SECTION  VIII 


LOAD  EXCEEDANCE  CURVES 


The  means  for  establishing  overall  load  exceedance  curves 
based  on  mission  considerations  are  outlined  in  this  section. 
These  exceedance  curves  may  be  desired  for  possible  use  in 
structural  fatigue  studies  of  the  airplane,  or  they  can  be  of 
value  in  comparing  how  the  load  exceedance  curve  of  one  mission 
compares  relatively  with  that  of  another.  Three  cases  differing 
slightly  in  detailed  makeup,  but  all  of  which  should  yield 
essentially  the  same  load  exceedance  curves,  are  described. 
Specifically,  the  cases  have  been  tailored  so  that  the  exceedance 
curves  would  pass  through  control  points  which  are  indicated  by 
equation  (25);  these  control  points  are  given  in  the  following 
listing: 


1000  ft 

P 

vfps 

X 

A 

Lxio1 

N0 

0 

.0005 

6 

54 

6.17 

5 

.0005 

6.6 

59.3 

6.17 

10 

.0005 

7.2 

64.8 

6.17 

20 

.0005 

8.4 

75.6 

6.17 

30 

.0005 

9.6 

86.4 

6.17 

35 

.0005 

10.2 

91.8 

6.17 

4o 

.000257 

8.7 

78.3 

3.17 

50 

.0000714 

5.8 

52.2 

.882 

Basically,  a  constant  value  of  =  9  was  chosen  for  all  the 
control  points.  °x 

Case  (a).-  Suggested  generalized  load  exceedance  curves  for 
this  case  are  given  by  the  following  equation 


_al  To~  ~a2  To ' 

=  .993  e  "  +  .007  e  w  (42) 


where 


-  V? 


a1  =  V.993  +  .007a 

a  =  h 

2  a 


This  equation  represents  a  modified  and  simplified  form  of  Cases 
k  and  £  presented  in  reference  2,  It  is  seen  to  be  a  function 
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of  just  three  environmental  parameters,  a  proportion  of  time  P  , 
an  intensity  a  ,  and  a  shape  parameter  a  ,  each  a  function  of 
altitude  only.  ^Recommended  values  of  these  three  parameters  are 
given  in  figure  29.  Basic  nondimen sional  exceedance  curves  de¬ 
rived  from  equation  (42)  using  the  values  of  a  given  in  figure 
29  are  shown  i.-  figure  34.  The  values  of  a  were  established  so 
that,  using  the  P  and  a  values  shown  in  figure  29,  the  ex¬ 
ceedance  curves  would  passwthrough  the  control  points  listed  at 
the  beginning  of  this  section.  The  use  also  of  the  values  of  P 
and  ct  given  in  figure  33  allows  the  curves  of  figure  30  to  be 
converted  into  the  convenient  working  form  shown  in  figure  31. 
Steps  for  establishing  the  expected  load  exceedance  curve  are  as 
follows : 


1.  Lay  out  the  mission,  as  in  figure  24,  and  establish  the 
A  and  N  value  for  each  segment  through  use  of 
equations°(15  and  (16). 

2.  Determine  the  composite  exceedance  curve  through 
either  of  the  following  equivalent  expressions 


n 


-  piTr 


P2T2N0, 


f  ( .  x  ^ 

>^W 


/ 


(*3) 


"  -  Vo^)  +  Vofffy  +  ■■■  ( 

where  f.  refers  to  the  curves  of  figure  30 
an  altitude  appropriate  to  the  segment  i  ,  while 
g*  refers  to  figure  31.  Illustrative  load  exceed¬ 
ance  curves  that  result  from  this  procedure,  and  the 
means  for  making  a  design  check,  if  desired,  are 
shown  in  figure  32.  For  example,  in  reference  to 
this  figure,  if  mission  A  is  taken  as  the  design 
mission,  and  if  the  design  requirement  is  stipulated 
that  the  number  of  positive  limit  load  exceedances 
must  not  exceed  5  in  30,000  hours  of  flight,  then 
the  curve  for  mission  A  for  30,000  hours  must  be 
such  that  It  falls  on  or  below  the  design  limit 
point  designated  by  P  on  the  figure. 


Case  (b).-  Although  not  essential,  generalized  exceedance 
curves  may  be  constructed  through  use  of  the  concept  of  dividing 
turbulence  encounter  into  two  parts,  one  associated  with  mild  to 
Intermittent  intensity  levels,  a  second  with  severe  levels. 
Because  of  the  popularity  of  this  concept  in  the  past,  curves 
based  on  this  concept  are  also  included  here.  The  specific 
equation  for  this  case  is 


+  P2e 


(45) 


74 


The  P  and  a  values  that  are  recommended  are  given  in  figures 
33  and  34.  From  these  values,  the  generalized  curves  of  figures 
35  and  36  follow,  where  for  figure  35 


P  "  P1  +  P2 

2  _  P12  ^2 

aw  "  P  °1  +  P 


Or 


(46) 

(*7) 


Again,  the  generalized  exceedance  curves  are  such  that  they  pass 
through  the  control  points  listed  at  the  beginning  of  this  section. 
The  values  of  P  and  a  as  found  from  equations  (46)  and  (47) 
are  very  similar  to  the  values  of  P  and  0  given  in  figure  29. 
Thus,  the  P  and  c  values  of  Case,  (a),  figure  29,  serve  also 

for  use  with  the  generalized  curves  for  this  case  given  in  figure 
35.  Note  that  the  curves  of  figures  35  and  36  agree  substan¬ 
tially  with  the  curves  of  figures  30  and  31.  Expected  load 
exceedance  curves  for  a  specified  mission  are  found  for  this  case 
in  the  same  manner  as  outlined  in  Case  (a).  It  is  to  be  noted 
that  the  cr2  and  Pp  values  of  this  case  are  simply  the  a  and 
P  values  given  in  Section  VII  on  design  procedures. 

Case  (c).-  This  case  is  included  to  show  that  certain  prop- 
erties  of  generalized  exceedance  curves  are  still  not  known,  and 
further  that  there  is  a  sensitive  relationship  between  the  shape 
of  the  curves  and  turbulence  intensity  values.  Each  portion  of 
the  curves  of  figure  35  is  derived  on  the  basis  that  the  r.m.s. 
value,?  or  turbulence  intensity  are  distributed  according  to  an 
equation  of  the  form 


p(°>  -  t  yf 


1  a 

2T 


This  equation  leads  to  an  exceedance  curve  of  the  form  e  x  , 
illustrated  by  curve  a  in  figure  37;  the  curve  was  readily 
accepted  and  became  popular  mainly  because  of  the  ease  of  handling 
the  simple  exponential  form.  There  is  evidence,  however,  that  the 
p(a)  distribution  may  be  more  in  the  form  of  a  Rayleigh- type  dis¬ 
tribution  given  by 


P (o)  =  2 


a 


This  distribution  leads  to  exceedance  curve  b  of  figure  37, 
see  Case  h,  figure  1,  of  reference  2;  the  equation  for  the  curve 
is 


N 


K 


(48) 
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t 

It  is  noted  that  experimental  data  often  exhibit  the  convex  shape 
noted  in  the  b  curve  for  low  values  of  the  abscissa,  thv.v  add¬ 
ing  substance  to  the  thought  that  the  p(a)  distribution  leading 
to  b  might  be  more  realistic.  With  the  use  of  equation  (48) 
an  equation  analogous  to  equation  (42)  may  be  written  as  follows 


W-Q  -  *1  -O'  Kl(^  al  +  a2  1%  Kl(^ 

where  also 

a1  ="^993  +  007a2 


Exceedance  curves  given  by  equation  (49)  and  the  values  of  a 
from  figure  29  are  shown  in  figure  38.  These  curves  may  be 
used  in  the  same  manner  as  for  Gases  (a)  and  (b)  but  they  call 
for  the  use  of  different  a  values,  specifically  those  re¬ 
commended  ir.  figure  39.  The  P  values  remain  the  same  as  given 
in  figure  29,  With  the  0  value  of  figure  39  and  the  P 
values  of  figure  29,  curves  quite  similar  to  those  shown  in  fig- 
ares  31  and  36  result.  Thus,  the  difference  in  shape  of  the 
curves  of  figure  38  relative  to  the  curves  of  figures  30  and 
35  is  compensated  for,  to  a  large  extent,  by  the  increased  a 
values. 

Stationary  aspects.-  Some  comments  relative  to  the  stationary 
properties  of  gust encounter  are  given  to  close  this  section.  The 
treatment  given  In  reference  2  indicated  that  generalized  exceed¬ 
ance  curves  could  be  derived  whether  the  r.m.s.  value  of  the  gusts 
was  stationary  or  not  (see  discussion  associated  with  equation  (11) 
in  reference  2) .  It  may  be  of  interest  to  demonstrate  the  con¬ 
struction  of  a  gust  response  model  that  is  implied  by  the  gener¬ 
alized  exceedance  curves,  which  shows  a  stationary  property  with 
respect  to  Nn  .  Consider  gust  encounter  in  terms  of  the  dis¬ 
crete  patch  concept  as  depicted  by  the  top  sketch  in  figure  40$ 
in  general,  the  r.m.s.  and  Nq  values  are  different  in  each 
patch.  (The  t*  values  are  assumed  to  correspond  to  the  P-jT.^ 
values  used  elsewhere  in  this  report.)  Various  relations  associa¬ 
ted  with  the  response  shown  are 


*2 

m 


t  =  t^  +  tg  +  t^ 


tN~  “  tlN01  +  t2N02  +  t3N03 


0 

2 

7x 


=  K 


.  2  ,  2  ,  ,  2 

Vi +  V2  +  Vv 


(50) 
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(51) 


n  -  Vo/lfe)  +  +  t3N03f^%) 

If  the  equation  for  n  is  multiplied  by  x  and  integrated,  in 
accordance  with  equation  (14)  of  reference  2,  then  the  following 
result  for  a  is  indicated 

■X 


^tiNo1ai 


fc2N0 


fc3N0, 


O' 


(52) 


This  result  appears  to  contradict  the  value  give*.,  by  equation 
(50).  Equation  (52) ,  however,  has  a  meaningful  interpretation. 
Assume  the  time  scale  of  each  patch  is  expanded  or  contracted  so 
thajt  the  Kg  value  for  each  patch  is  the  same,  specifically  N0  • 
The  change  is  easily  made  by  the  transformation 

*  \ 

*1  ■  ‘i 


The  transformed  system  Is  depicted  by  the  bottom  sketch  of  figure 
40.  It  is  noted  that  the  transformation  does  not  change  the 
r.m.s.  value  of  each  patch,  nor  does  it  change  the  peak  count  of 
the  record.  The  mean  square  value  of  the  new  system,  however,  is 
changed  and  Is  precisely  that  given  by  equation  (52).  It  Is 
thus  seen  that  the  construction  of  exceedance  curves  in  the  format 
given  by  equations  (43)  or  (51)  and  use  of  a  mean  square  defin¬ 

ition  given  by  equation  (14)  of  reference  2  effectively  converts 
the  response  record  into  one  which  is  stationary  with  respect  to 
N0  .  Generally,  there  will  be  little  difference  between  the 
r.m.s.  values  given  by  equations  (50)  and  (52).  Equation  (50) 
rather  than  equation  (52)  was  therefore  adopted  as  the  basis 
for  establishing  the  composite  r.m.s.  values  of  Section  VTI 
(equation  (38)),  since  it  is  somewhat  easier  to  evaluate. 
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SECTION  IX 


SPECTRAL  VS.  DISCRETE-GUST  DESIGN 


The  concern  of  some  is  that  while  spectral  approaches  to  the 
gust  problem  may  be  appropriate  in  many  respects,  they  may  not  be 
so  in  certain  respects  relative  to  the  discrete-gust  design  method. 
One  concern  is  that  the  spectral  approaches  are  more  involved 
computationally;  another  is  that  the  spectral  approach  may  not  in¬ 
clude  the  consideration  of  the  severe  isolated-type  gust  that 
seemingly  is  encountered  on  occasions.  The  purpose  of  this  sec¬ 
tion  is  to  show  that  the  spectral  approaches  of  this  manual 
effectively  obviate  both  of  these  concerns.  Indeed,  particular 
note  should  be  taken  of  the  fact  that  the  preliminary  procedure 
recommended  herein  is  quite  equivalent  to  the  discrete-gust  method; 
it  is  just  as  simple,  has  the  same  basic  steps,  but  at  the  same 
time  includes  in  better  fashion  inherent  dynamic -behavior  effects 
of  the  airplane. 

It  should  be  noted  that  in  the  discrete-gust  approach,  the 
choice  of  the  1-cosine  gust  has  really  no  sound  basis.  A  sine 
gust,  a  triangular  gust,  or  even  a  rectangular  gust  would  serve 
as  well.  Thus,  the  representation  of  gust  encounter  by  a  single 
type  gust  is  quite  an  arbitrary  matter.  It  would  seem  that  the 
question  of  Mgust  shape"  in  the  discrete-gust  approach  poses  just 
as  great  a  problem  as  is  involved  in  deciding  whether  the  gusts 
should  be  random  or  of  discrete  shape. 

For  most  gust  encounter  there  is  little  doubt  that  the  pro¬ 
cess  is  continuous  and  random.  Spectral  analyses  are  thus  appro¬ 
priate.  It  may  be  asked,  however,  whether  the  spectral  approach 
includes  the  effects  of  large  isolated  gusts,  assuming  they  do 
exist.  The  intent  of  the  analysis  leading  to  equation  (11), 
reference  2  is  to  take  into  account  such  possibilities.  Thus,  if 
the  function  g  used  in  this  equation  is  chosen  properly  and 
contains  the  Information  relative  to  the  large  gusts,  then  these 
gusts  should  be  reflected  in  the  peak  count  curve.  Further,  gust 
intensity  values  out  to  infinity  are  usually  considered  (for 
mathematical  convenience)  in  the  derivation  of  the  exceedance 
curves.  The  tail  portion  of  the  exceedance  curves,  thus,  intro¬ 
duces  extreme  intensities,  even  though  they  do  not  exist  in  prac¬ 
tice.  The  effect  on  the  exceedance  curves  of  placing  a  finite 
upper  limit  on  the  gust  intensity  that  can  be  encountered  is  shown 
in  the  section  on  ultimate  load  (Section  X) . 

Some  of  the  thoughts  of  the  preceding  paragraph  may  be  illus¬ 
trated  more  specifically  by  the  following  analysis  which  is  intend¬ 
ed  to  show  that  the  preliminary  spectral  approach  given  herein  is 
quite  comparable  to  the  discrete  approach.  The  acceleration  given 
by  ohe  discrete  approach  is 


Multiply  equation  (22)  by  n  to  obtain 

n0An  =  K0H°w  =  Arjaw 

Consider  now  that  rja^  equals  An  ,  then  these  two  equations 


v  Uj  .  =  no' 
de  1  w 

SoXiXSS  ^elocity  value  no  is  thus  established  that  may  be 

does1^^  in  wh%Same  ?°le  in  the  spectral  approach  as 
cation  of  ^  dj;crete-pst  approach;  that  is,  the  multipli¬ 
cation  of  qaw  by  the  spectral  value  of  A  yields  the  same 
acceleration  as  does  the  discrete-gust  approach  Val np «  nf*  - 
as  obtained  from  this  relation  andSthe  K^ataei  o?  ftarf  16 
are  shown  in  figure  41.  using  Ude  *  50%s  .  This  010^1^^ 

lePvel  JSthS  lalllTof  ,  Wj  g.  ?'ua§a 

are^pecified'in  Sis“f  M  IT™  ^aSs  a°?  ^  0 

oreis3t?ur  while  Pua°e  °n  00mparable  terms  wlth  ude  since 
S'  la  'S’  hJ-le  is  an  equivalent  gust  velocity.  Based  on 

fb®  V;  3  value  given  in  figure  29,  the  choice  of  noT-54 
indicates  a  value  of  q  -  18  ;  with  reference  to  the  oIT"  6 

desivn  at  ii6ure  21,  the  value  which  is  critical  for 

num^sMst^  ift^Vt^T0^  ?n=lad®“inn«-eacomiiSltSe:rnt 

numDers  listed  in  the  introductory  remarks  of  Section  J. 
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Figure  41.  Spectral  "Velocities"  Equivalent  to  the 
Discrete-Gust  Design  Velocity 


SECTION  X 


VIEWS  ON  ULTIMATE  LOAD  DESIGN 


Gust  design  loads  associated  with  limit  load  levels  are 
recommended  herein;  gust  design  in  terms  of  ultimate  loads  is  not. 
The  main  reason  is  because  of  the  large  uncertainty  of  the  load 
exceedance  values  in  the  region  of  ultimate  load  levels.  With 
the  presently  available  exceedance  curves,  there  is  still  some 
uncertainty  in  the  range  of  limit  load  values  (x/ax  =  18  to  30), 
and  the  uncertainty  becomes  ever  increasing  for  load  levels  that 
extend  out  to  the  ultimate  load  range  (x/ax  =  27  to  45  ,  on  the 
assumption  that  ultimate  strength  levels  are  1.5  limit  strength). 
In  fact,  data  on  load  exceedance  in  the  ultimate  load  range  are 
virtually  nonexistent. 


The  marked  sensitivity  of  the  load  exceedance  curves  in  the 
region  of  large  load  levels  can  be  demonstrated  by  studying  the 
influence  of  various  finite  limits  on  representative  joint- 
probability  distribution  functions  for  gust  encounter.  Specific¬ 
ally,  a  study  of  the  effect  of  using  a  finite  upper  limit  in 
equation  (11a)  of  reference  2  is  pertinent.  Normally,  exceedance 
curves,  as  established  in  reference  2,  are  derived  on  the  assump¬ 
tion  that  the  values  of  x  and  x  in  the  joint  distribution 
function  are  unbounded;  that  is,  they  extend  to  infinity.  This 
assumption,  of  course,  is  not  realistic.  Nature,  itself,  seems 
to  have  a  self-limiting  process  on  the  maximum  gust  velocity  that 
can  develop.  If  this  value  is,  say,  300  fps  ,  and  if  the  overall 
r.m.s.  value  of  turbulence  encounter  is  3  fps  ,  then  it  would 
appear  that  x/ox  cannot  ever  be  larger  than  100  .  Other 

limiting  processes  also  enter.  For  example,  a  gust  velocity  of 
100  fps  may  cause  the  airplane  to  stall.  Thus,  the  maximum  load- 
producing  velocities  may  be  limited  to  100  ,  indicating  a  maximum 
value  of  x/ax  s  30  . 

To  study  how  the  tails  (or  absence  of  tails)  of  the  joint- 
distribution  curve  affect  the  exceedance  curves,  equation  (ll) 
of  reference  2  was  used  in  conjunction  with  some  of  the  Joint 
distribution  curves  also  given  in  that  reference,  with  the  restric¬ 
tion  that  the  equation  applied  only  within  the  region  defined  by 


Outside  of  this  limit,  p  was  taken  as  zero.  It  was  found  that 
the  exceedance  curves  so  obtained  could  be  approximated  with  good 
accuracy  by  simply  subtracting  a  constant  value  from  the  exceedance 
curves  that  were  derived  with  an  upper  limit  of  infinity.  The  con¬ 
stant  is  chosen  simply  as  the  value  of  the  unbounded  exceedance 
curve  at  the  abscissa  equal  to  the  upper  limit  bound  of  x/cx  =  R; 
the  value  of  the  bounded  exceedance  curve  thus  approaches  zero  as 
the  limit  R  is  approached.  Some  results  associated  with  the 
exceedance  curves  of  Case  (-c-)>  figure  30,  are  shown  in  figure 

42.  in  the  region  of  large  x  ,  it  is  seen  that  the  exceedance 
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values  are  very  sensitive  to  the  cutoff  value  R  .  It  may  be 
observed  that  the  curves  for  finite  R  are  more  realistic  than 
the  curves  for  R  =  oo  ,  since  the  concept  of  a  joint-probability 
distribution,  confined  to  within  a  finite  region,  is  probably  more 
realistic  of  nature.  There  is  no  way  to  ascertain,  however,  which 
of  the  curves,  if  any,  most  closely  represent  realistic  behavior. 
This  large  uncertainty,  then,  is  one  of  the  reasons  why  it  is  felt 
appropriate  to  avoid  ultimate  load  design  concepts. 
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Figure  42.  Effect  of  Bounded  Joint  Distribution 
on  Exceedance  Curves 
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Previous  sections  of  this  manual  have  dealt  mainly  with 
response  quantities  that  can  be  treated  in  an  independent  or  sep¬ 
arate  manner.  Often,  however,  a  structural  component  is  designed 
by  loads  or  stresses  which  act  in  combination.  The  manner  of 
handling  such  situations  is  described  in  this  section.  There  are 
two  basic  combined  loads  problems;  each  will  be  treated  separately. 

Stresses  due  to  superposition.-  Reference  4  indicates  the 
means  for  handling  this  case.  Insummary,  the  stress  at  a  crit¬ 
ical  point  in  the  structure  may  be  due  to  several  loads  that  act 
simultaneously,  for  example 

s  =  ax  M(t)  +  a2V(t)  +  a3T(t) 

where  M  ,  V  ,  and  T  refer  to  moment,  shear,  and  torque,  respec¬ 
tively.  The  proper  way  to  handle  this  case  is  to  establish  the 
frequency  response  function  for  s  ,  thus 

H8(a>)  =  ax  Hm(o>)  +  a2Hv(o>)  +  a^H^o) 

This  frequency  response  function  is  used  to  establish  the  power 
spectrum 

<p s H  -  Ihs(o>)12  $w(o>) 

Treatment  thereafter  is  in  the  same  manner  as  for  an  individual 
stress.  By  this  approach,  all  phase  relationships  between  the 
moment,  shear,  and  torque  will  be  taken  into  account  automatically. 

Combined  stresses  by  interaction.-  A  general  treatment  of 
combined  stresses  by  Interaction  is  given  in  reference  5.  Only  a 
summary  of  some  of  the  simplified  aspects  of  the  treatment  will  be 
given  here.  If  a  design  check  fails  to  meet  the  tests  that  are 
outlined  herein,  then  recourse  should  be  made  to  the  more  detailed 
procedures  given  in  reference  5. 

Consider  that  a  structural  element  is  under  the  action  of 
two  stresses,  such  as  an  axial  stress  and  a  shear  stress,  and  that 
the  integrity  of  the  element  Is  governed  by  an  interaction  curve 
in  terms  of  these  two  stresses.  The  problem  is  to  find  the 
expected  number  of  crossings  of  the  interaction  border.  Figure 
43  illustrates  the  case. 

Steps  for  making  a  simplified  solution  based  on  using  an 
approximating  elliptical  interaction  curve  are  as  follows: 

1.  As  in  the  preliminary  approach,  decide  on  the 
critical  flight  conditions. 

2.  By  equations  (8)  to  (15)  compute  as  ,  oT  , 

0s  »  aT  >  P  >  an<3  v  >  note,  for  convenience  choose 
ow  =  1 
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3-  By  the  sequential  graphical  construction  shown  in 
figure  4 4,  establish: 

a)  cfsm  *  and  from  as  >  i  a™3  p; 

since  aw  =  1  was  used  in  step  2,  then  As  =  crs 

and  AT  =  oT  (top  sketch), 
m  Tm 

b)  °hQ  *  and  vo  f’rom  as  >  aT  >  v  * 

and  0O  of  step  3(a)  (middle  sketch). 


e)  sa  from  os  ,  aT  ,  0£  ,  at  and 

v0  of  steps  3(a)  ana  3(b)  (bottom  sketch). 


4.  Refer  to  figure  45.  Through  origin  of  s  ,  r  axes , 
which  is  located  at  the  1-g  load  point,  draw  sm  and 
Tm  axes;  the  direction  of  rotation  from  s  to  sm 
is  the  same  as  the  direction  of  rotation  of  diameter 
A-B  in  figure  44  to  the  asm  and  aTm  values.  Estab¬ 
lish  ellipses  on  the  sm  ,  Tm  axes  according  to 
the  equation 


for  various  values  of  z  . 


5.  Choose  the  largest  ellipse  which  Just  begins  to 
exceed  the  boundaries  set  by  the  interaction  curve 
and  note  the  value  for  z  for  this  ellipse. 


6 .  Prom  os  ,  oT 


“m  ‘m 

3,  determine  N0  by  the  equation 

„  1  sa 

N°  -5F  oTS— 

m  m 


of  step  1,  and  from  s  of  step 

cl 


7.  With  the  value  of  z  established  in  step  5,  the 
value  of  N0  of  step  6,  the  and  P  values 

of  figures  21  and  22,  and  a  chosen  time  of  flight 
T  ,  establish  the  number  of  expected  exceedances  by 
the  equation 


z 

n  =  1.596  PTN0  /£.  p 

This  equation  represents  the  approximation  of 
equation  (38),  n  =  1/2  ,  of  reference  5,  for 
large  z 
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8.  Judge  whether  the  value  for  n  obtained  in  step 

7  is  acceptable;  the  acceptable  value  of  n  depends, 
of  course,  on  the  nature  of  the  interaction  curve  - 
whether  it  is  a  strength-type  boundary,  a  buckling 
boundary,  or  some  other  form.  If  n  is  judged  too 
large,  proceed  to  step  9- 

9.  If  n  of  step  7  is  judged  to  be  too  large,  then 
proceed  as  follows: 

a)  In  terms  of  the  elliptic  boundary  approach 
described  in  the  previous  steps,  re-evaluate  n 
by  means  of  a  mission  segment  consideration# 
similar  to  that  given  in  the  section  on  detailed 
design  considerations;  thg,t  is,  use  the  equation 
in  step  7  in  a  way  analogous  to  that  indicated  by 
equation  (4l) . 

b)  If  n  as  found  in  (a)  is  still  judged  too  large, 
then  recourse  should  be  made  to  the  detailed 
method  of  treatment  that  is  outlined  in  reference 
5,  using  the  given  interaction  curve,  in  segmented 
fashion.  The  use  of  the  elliptic  curve  approxi¬ 
mation  yields  an  upper  bound;  in  general,  the  use 
of  the. actual  interaction  curve  will  lead  to 
smaller  values  for  n  . 
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Figure  4^.  Approximating  Elliptical  Interaction  Boundary 
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SECTION  XII 


^  NONUNIFORM  SPANWISE  GUSTS 


This  section;  outlines  the  ;means  for  treating  the  gust 
encounter  case  wherein  the  gusts  are  considered  random  in  the 
spanwise  directionsf  as  wellas4n  the  direction  of  flight,  as 
depicted  in  the  fbllowiqg^sketch.  < 


This  case  may  be  of  significance  for  large  flexible  aircraft  where 
the  span  may  become  ah  appreciable  fraction  of  the  scale  of  turb¬ 
ulence.  In  these  circumstances,  it  is  possible  for  some  of  the 
modes,  particularly  the  antisymmetrical  and  higher  frequency  modes, 
to  be  excited  to  a  much  greater  extent  than  they  would  be  in  the 
case  of  uniform  spanwise  gusts.  The  overall  response  for  the 
completely  random  gust  environment  may  thus  be  larger  than  for  the 
environment  which  is  considered  random  only  in  the  flight  direc¬ 
tion. 


Detailed  consideration  of  this  case  is  given  in  reference  29. 
Essentially,  the  case  is  that  of  a  system  which  is  being  excited 
by  multiple  random  inputs  rather  than  by  a  single  one.  As  such, 
the  cross  spectra  between  the  individual  inputs  must  be  considered, 
along  with  the  spectral  values  of  each  input.  The  cross  spectra 
which  are  of  concern  here,  and  which  are  related  to  the  basic  gust 
input  spectrum  of  this  manual,  equation  (17)  >  have  been  estab¬ 
lished  by  the  principal  investigator,  reference  30  .  Some  of  the 
basic  derived  results  are  shown  in  figure  46.  The  means  for 
using  these  results  in  gust  response  calculations  is  quite  straigit- 
forward,  involving  mainly  an  increase  in  computational  effort. 
Suggested  procedures  based  on  the  work  of  reference  29  are  as 
follows: 


1.  Divide  the  wing  into  equally  spaced  spanwise  seg¬ 
ments;  from  6  to  10  segments  over  the  span  should 
be  adequate  for  most  cases.  Suppose  the  length  in 
the  span  direction  for  each  segment  is  designated 
by  e  . 


For  each  segment,  derive  the  frequency  response 
functions  for  the  response  variable  of  concern, 
considering  that  the  sinusoidal  gust  encounter 
acts  only  over  the  segment  being  treated;  this  is 
similar  to  a  strip  theory-type  of  treatment,  where* 
in  each  segment  is  handled  separately.  Call  the 
frequency  response  function  ,  where  m  desig¬ 
nates  the  number  of  the  segment. 

With  the  Hm  functions  of  step  2,  and  the  cross 
spectra  of  figure  46,  establish  the  spectrum 
for  response  according  to  the  equation 

0X  =  0oUHil2  +  1h2'2  +  1h3|2  +  •••] 


+  0Q1  2  Re (H-jHg  +  H2H3  +  ...) 
+  <t>02  2  Re +  H2H4  +  ...) 


+  •  • 


+  0Qm  2  Re  (H  A  +  ¥m+l  +  *  •  •  ) 

where  0  is  the  spectrum  for  £  *  0  (the  same 
u  as  equation  (17) 

0^  is  the  cross  spectrum  for  I*  =  q/L  =  me/L 
and  where  Re  denotes  the  real  part. 

4.  With  0X  of  step  3,  establish  A  and  NQ  as 
with  any  output  spectrum. 

5.  Use  the  A  and  N  values  of  step  4  in  the  same 
manner  as  used  in  other  sections  of  the  manual. 


001  .01  .1  I.  (0  100 


Fisure  46,  Cross-Spectra  for  Treatment  of  Nonuniform 
Spanwise  Gusts 


APPENDIX  A 


STANDARD  DATA  AND  CONVERSION  CHARTS 


This  appendix  includes  data  which  may  be  useful  in  studying 
various  gust  response  problems.  Included  are  some  data  on  a 
standard  atmosphere,  some  graphs  pertaining  to  common  probability 
density  and  probability  functions,  and  several  conversion  graphs 
related  to  frequency  and  wavelength.  Various  definitions  involved 
in  the  conversion  charts  are  the  following: 


iL  -  10  =  V  _  Vk 
2tt  2ir  T.  ttc 


to  --  2rrf  =  Vft  =  =  — 

A.  C 

n  27Tf  to  2tt  2k 

“  =  ^r  =  v  =  x  =  T 

.  TTfC  COC  ftc  TTC 

k  =  —  =  2V  =  T  =  'T 

to  -  2irLi '  Lto  _  2irL  _  2Lk 

- -  ^ - ~ 

>  _  V  _  2ttV  27T  ire 
*  “  7  to”  =  TT  “  TT 


where 


TABLE  I  -  SOME  STANDARD  ATMOSPHERIC  PROPERTIES 


h 

1000  ft 

T 

°R 

p  X  10^ 
slugs/cu  ft 

p/p0 

c 

ft/sec 

0 

518.7 

2.377 

1.000 

1117 

5 

500.9 

2.048 

.8617 

1098 

10 

483.0 

1.755 

.7385 

1078 

15 

465-2 

1.496 

.6292 

1058 

20 

447.4 

1.266 

.5328 

1037 

25 

429-5 

I.O65 

.4481 

1016 

30 

411.7 

.8893 

•  3741 

995.1 

35 

393.9 

•7365 

.3099 

973.3 

36.089  390.0 

.7061 

.2971 

968.5 

40 

390.0 

.5851 

.2462 

968.5 

45 

390.0 

.4601 

.1936 

968.5 

50 

390.0 

.3618 

.1522 

968.5 

55 

390.0 

.2845 

.1197 

968.5 

60 

390.0 

.2238 

.09414 

968.5 

65 

390.0 

.1760 

.07403 

968.5 

Note:  here  T  represents  temperature 
and  c  speed  of  sound. 
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Figure  A-2.  Probability  Values  for  Curves  of  Fig.  A-l 


f  ,cps 

Figure  A-3.  Nomograph  Relating  Q  and  f 
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i,-^KTiVr^StyiKJ'/iAA' 


inn  nnn 


0001  .001  .01 


&,  radians/  ft. 


Figure  A-7.  Nomograph  Relating  X  arid  ft 
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APPENDIX  B 


DETERMINATION  OF  FREQUENCY  RESPONSE  FUNCTIONS 


Since  the  frequency  response  function  plays  such  a  central 
role  in  gust  reponse  analysis  by  power  spectral  methods,  a  review 
is  given  in  this  appendix  of  some  of  the  ways  for  establishing  this 
function.  Throughout  the  presentation,  sinusoidal  motion  is 
implied . 

Structural  aspects.-  Two  prime  means  are  normally  used  to  treat 
the  aircraft  structure,  one  a  modal  approach,  the  other  a  lumped- 
mass  approach . 

Modal  approach: 

The  response  or  deformation  of  the  structure  is  assumed  to  be 
given  in  terms  of  the  natural  modes  of  the  structure,  such  as 


Z  =  a1Z1  +  a2Z2  +  a^Z^  +  ...  (B-l) 

where  the  Z's  are  the  modes,  and  the  a's  are  coefficients  to  be 
established.  A  Lagrangian  dynamical  treatment  leads  to  the  equa¬ 
tion 


HnSn  +  Vfcn  ’  Jf  * n  as 


(B-2) 


where  M^  is  the  generalized  mass  of  the  Z„  mode  ,  Mi  =  //mZ2ds. 

is  the  natural  frequency  of  the  mode  Z«  ,  and  71  JJ  n 

p  is  the  applied  loading  over  the  surface  S  .  Structural  damp¬ 
ing  is  included  by  adding  a  term  27nMn<Mn  to  the  left  side,  or 
by  modifying  the  second  term  to  read  Mn  o$(H-ign)an  . 

Lumped-mass  approach; 

The  mass  of  the  structure  is  represented  by  "equivalent"  con¬ 
centrated  masses;  the  stiffness  is  represented  by  discrete  springs, 
or  the  concept  of  influence  coefficients  is  introduced.  The 
response  equation  takes  the  form  of  a  matrix  equation 


[DjJZgi  =  o>2['m-]jZs|  +  |  Ps  | 


(B-3) 


where  Ifc  represents  applied  concentrated  loads.  Structural  damp¬ 
ing  is  often  introduced  by  multiplying  the  D  matrix  by  a  term 
(1+ig)  • 

Aerodynamics,.  The  aerodynamics  are  handled  usually  by  a  strip 
theory  approach  or  by  lifting  surface  theory;  only  the  lifting 
surface  theory  method  will  be  discussed  here,  and  only  for  sub¬ 
sonic  flow.  Two  forms  of  lifting  surface  theory  have  been  used, 
one  which  uses  a  kernel  function  and  assumed  pressure  modal  shapes, 

the  other  which  uses  a  modified  kernel  form  and  concentrated  aero¬ 
dynamic  loads. 
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The  basic  kernel  approach; 


The  downwash  w  that  is  produced  by  ah  aerodynamic  pressure 
distribution  p  is  given  by  ihe  equation 

w(x,y)  =  jj  p(£,T})K(x-S,y-Ti)dx  dy  (B-4) 


where  K  is  a  kernel  function  which,  analogous  to  an  influence 
coefficient,  gives  the  downwash  at  point  x,y  due  to  a  unit 
pressure  dipole  at  £,t]  j  it  is  a  function  of  reduced  frequency 
k  =  oxj/2V  ,  and  Mach  number  M  ;  see,  for  example,  reference  31  • 

For  a  surface  which  is  oscillating  with  a  deflection  form 
Z  and  for  the  condition  of  sinusoidal  gust  encounter,  the  down- 
wash  is  given  by 


w  =  V 


dZ 

3x 


-i 

+  io£  -  e 


cox 

V 


(B-5) 


where  x  is  the  location  of  the  downwash  point  relative  to  some 
convenient  reference  point  such  as  the  leading  edge  of  the  root 
chord  of  the- wing.  The  problem  is  to  combine  equations  (B-4)  and 
(B-5)  and  to  solve  for  .the  pressure  p  due  to  the  motion  Z  and 
the  gust  input  e"i\aacA7  , 

Modified  kernel  approach: 

In  the  modified  kernel  approach  the  aerodynamic  loading  is 
assumed  to  be  replaced  by  equivalent  concentrated  loads.  In  terns 
of  these  loads,  equation  (B-4)  may  be  converted,  see  reference  3, 
to  the  matrix  form 


iw|  =  [K]|P|  (B-6) 


where  K  is  a  square  matrix  composed  of  terms  of  a  modified  ker¬ 
nel  function,  developed  in  reference  3*  Reference  32  presents  a 
related  treatment  which  also  leads  to  an  equation  similar  to 
equation  (B-6).  The  downwash  due  to  the  lifting  surface  and  the 
gust  is  the  same  as  that  given  by  equation  (B-4);  for  this  case, 
solution  for  P  can  be  made  directly  because  of  the  fact  that 
the  iC  matrix  of  equation  (B-6)  can  be  inverted. 

Frequency  response  function.  -  Four  means  are  outlined  now  for 
establishing  the  frequency  response  function.  One  uses  the  ker¬ 
nel  function  approach)  the  other  three  are  in  terms  of  the 
modified  kernel. 
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In  terms  of  the  kernel  function; 

Substitute  equation  (B-l)  into  equation  (B-5)  and  obtain 

imx 

w  =  a.^  +  a2z2  +  a^z^  +  ...  e  V  (B-7 ) 

/  d  \ 

where  zn  are  downwash  values  given  by  ( ' V  ^  +  io JZn  .  Assume 

p  =  b1p1  +  b2P2  +  b3P3  +  ...  (B-8) 

where  the  p's  are  assumed  pressure  mode  functions.  Substitute 
equation  (B-1^)  in  equation  (B-4),  evaluate,  and  obtain 


w  =  b1w1  +  b2w2  +  b3w3  (B-9) 

where  the  w^s  are  the  downwash  functions  associated  with  the 
Pns  .  At  a  number  of  points  on  the  wing,  equal  In  number  to  the 
number  of  b^s  that  are  chosen  in  equation  (b-8),  equate  the 
values  of  w  given  by  equations  (B-7)  and  (B-9);  this  leads  to 
a  matrix-type  equation  of  the  form 

^i(l)  w2 (l )  w3(l)  b^  z1(l)  Zg(l)  z3(l)  z^(l) 

w1(2)  w2(2)  w3(2)  b2  =  Z;l(2)  z2(2)  z3(2)  z4(2)  a2  c2 

+ 

w^(3)  ^2(3)  *3(3)  2^(3)  Zg(3)  z3(3 )  z4(3 )  a3  c3 

a4 

In  general,  the  number  of  an  values  is  different  from  the  number 
of  bn  values.  Invert  this  equation  to  obtain 

bl  “  Alal  +  ^2a2  +  A3a3  +  ^4a4  +  dl  ' 

b2  =  ....  I  (B-10) 

b-  =  . . . .  j 

j 

Substitute  equation  (B-8)  into  equation  (B-2),  evaluate  the  inte¬ 
grals,  then  make  use  of  equation  (B-10)  to  obtain  final  equations 
which  are  in  terms  of  the  unknown  a^s  only.  Solution  of  these 
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final  equations  for  aj^  allows  the  frequency  response  function 
for  various  quantities  of  interest  to  be  evaluated.  Generally, 
after  solution  is  made  for  the  an  values,  the  total  loading 
on  the  structure,  including  the  inertial  and  aerodynamic  loadings, 
is  used  to  compute  the  section  shear,  moment,  and  torque  loads. 
This  total  loading  is  given  by 


pt  =  mco2(a1Z1  +  a2Z2  +  ...)  +  p  (B-ll) 

where  p  is  given  by  equations  (B-8)  and  (B-10). 

By  the  modified  kernel  function: 

Scheme  1.-  This  scheme  makes  use  of  equation  (B— 3 ) j  in  general, 
the  load  and  deflection  points  used  for  the  dynamical  treatment 
of  the  system,  equation  (B-3),  may  not  be  the  same  as  those  used 
for  establishing  the  aerodynamic  lOcds,  equation  (B-6).  By  suit¬ 
able  interpolation  formulae,  however,  these  two  load  and  deflec¬ 
tion  systems  may  be  related;  suppose  that  the  transformations  from 
one  system  to  the  other  are  given  by 


where  T^ 
equations 


|P|  =  (T1]|Psl  (B-12) 

|Zal  -  tT2]|z|  (B-13) 

and  To  are,  in  general,  rectangular  matrices.  Combine 
(B-3),  (B-5),  (B-6),  (B-12)  and  (B-13)  to  obtain 

_i  oac 

K  TX(D  -  A)T2Z  “  (v  +  *»)  Z  -  e  V  (B-14) 


where,  for  simplicity  in  writing,  the  matrix  notation  has  been 
dropped.  Equation  (B-14),  in  turn,  may  be  written 


[k  t2(d  -  A)t2  -  (v  ^  +  i©)jz 


oac 

V 


or 


(B-15) 
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Inverted,  the  desired  result  for  the  frequency  response  function 
Z  is  found  to  be 


Prom  Z  ,  the  frequency  response  function  for  other  quantities  of 
interest,  such  as  acceleration  or  load  at  given  points,  may  be 
determined . 

The  equation  describing  flutter  of  the  system  follows  directly 
from  equation  (B-15),  being  simply  this  equation  with  the  right- 
hand  side  set  equal  to  zero. 

Scheme  2 .  -  If  equation  (B-2)  is  applied  to  each  of  the  modes  con¬ 
sidered  in  equation  (B-l),  and  if  the  loading  p  is  expressed  In 
terms  of  equivalent  concentrated  loads  at  various  grid  points, 
then  for  the  sinusoidal  case  of  an  the  following  matrix  equation 
can  be  written 


[(«£  -  “2Xj|an!  '  tZn3lpl  (B-17> 

where  -  u>2)MnJ  •5Ls  a  simPle  diagonal  matrix,  an  a  column 

matrix,  Zn  is  a  rectangular  matrix  built  up  from  row  matrices 
which  express  the  deflection  at  each  grid  point  for  each  mode  Zn, 
and  P  is  a  column  matrix  of  the  applied  concentrated  loads  at 
each  grid  point.  The  simple  matrix  form  of  the  generalized  forces 
on  the  right  side  of  equation  (B-l 7)  follows  from  equation  (B-2) 
because  of  the  Dirac  function  nature  of  each  concentrated  load. 
Now,  make  use  of  a  matrix  representation  of  equation  (B-l)  to 
express  the  deflection  at  each  of  the  chosen  grid  points,  namely 


|Z|  -  [Z^llaJ  (B-18) 

where  [z']  is  a  rectangular  matrix  in  which  each  column  repre¬ 
sents  thendef lection  at  the  various  grid  points  in  a  given  mode. 

Invert  equation  (B-6)  and  combine  with  equations  (B-5)  and 
(B-18)  to  obtain  the  loading  as 


|P|  - 


dx 


+  10) 


a„ 


n' 


-  tier1 


-i 


cox 

T1 


(B-19) 


The  substitution  of  this  equation  into  equation  (B-iy)  leads  to 
the  result 
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2^ 

or  Jm. 


'„j-  [ZJ[Srl[VS  +  HtZnl 


1_i  OJX 
e  v 


or 


i  -i  — ! 

[ElUJ  =  -  [Zn][K]"1|e  Vi 


(B-20) 


Solution  for  an  yields  the  frequency  response  functions  for  an  . 
With  the  a^s  established,  various  other  frequency  response  func¬ 
tions  such  as  acceleration  or  load  at  a  given  point  follow  readily. 

Scheme  3 .  -  In  tiis  version,  proceed  by  solving  equation  (Brl7) 
for  [anj  ,  or 


where  solution  follows  simply  because  the  square  matrix  on  the 
left-hand  side  of  equation  (B-17)  is  diagonal.  The  combination 
of  this  equation  and  equation  (B-l8)  yields 


CZn]|P| 


This  equation  and  equation  (B-5)  yield 


w 


(B-21 ) 


Combine  this  equation  with  equation  (B-6)  by  eliminating  w  and 
find  the  following  equation  in  terms  of  the  P  values  only 


COXi 

1  (B-22 ) 


or 
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[Hi  Ip!  «  - 


cox 

V 


Inverted,  the  frequency  response  values  for  P  are  found  directly 
as 


-i  a® 

lP|  =  -[H]'1  e  V  (B-23) 


It  is  to  he  noted  that  the  matrices  leading  to  H  are  all  very 
easy  to  establish  and,  thus,  this  version  appears  especially  attrac¬ 
tive;  notice  also  that  the  values  of  P  are  obtained  in  direct 
fashion. 

The  equations  for  flutter  follow  from  schemes  2  and  3  by 
simply  setting  the  right-hand  sides  of  equations  (B-20)  and  (B-23) 
equal  to  zero. 

Example  cases  of  frequency  response  functions  and  out put  spectra.- 
Several  frequency  response  functions,  as  obtained  in  the  study  of 
specific  aircraft,  are  given  here  to  illustrate  the  general  nature 
of_  the- functions.  Figures  B-l  and  B-2  show  representative  functions 
for  the  KC-135  airplane;  figure  B-3  for  the  B-58;  while  figure 
B-4  applies  to  the  B-57.  Two  representative  output  spectra,  ob¬ 
tained  by  multiplying  the  gust  spectrum  by  the  functions  of  fig¬ 
ures  B-l  and  B-2,  are  shown  in  figures  B-5  and  B-6. 

Some  comment  on  the  number  of  modal  functions  that  should  be 
included  in  any  response  treatment  is  considered  worthwhile.  In 
general,  it  is  not  possible  to  give  rigid  guidelines  on  how  many 
modes  should  be  included,  but  some  observations  found  in  practice 
give  certain  useful  hints.  For  the  larger  aircraft  it  is  usually 
found  that  the  range  of  frequencies  below  10  cps  is  of  concern; 
and,  more  often  than  not,  only  modes  below  5  cps  are  of  importance. 
The  mode  selection  is  tied  in  with  the  location  on  the  airplane 
where  the  response  is  being  considered.  For  bending  moment  near 
the  root,  for  example,  the  higher  modes  do  not  seem  to  contribute 
much.  On  the  other  hand,  the  higher  modes  are  often  found  to  con¬ 
tribute  substantially  when  response  near  the  outboard  regions  is 
being  determined.  Also,  when  main  loads  in  the  wing  are  being 
considered,  secondary  modes,  such  as  are  associated  with  pronounced 
pylon  lateral  motion,  generally  can  be  ignored.  In  computing  the 
loads  in  the  pylon,  however,  then  these  secondary- type  modes  become 
significant . 

The  main  reasons  why  the  higher  modes  (say  modes  above  3-4 
cps)  do  not  usually  contribute  significantly  to  the  response  is, 
of  couise,  due  to  the  rapid  fall-off  of  the  gust  input  spectrum 
in  the  high  frequency  regions.  Figures  B-2  and  B-6  illustrate  this 
effect  rather  well.  In  figure  B-2,  for  frequency  response,  it 
appears  that  four  higher  modes  are  of  concern  at  frequencies 
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of  3*3  >  4.2  ,  5 .8  and  6.6  cps  *  In  figure  B-6,  however,  it 
is  seen  that  the  two  higher  modes  have  almost  disappeared,  and 
that  the  response  contribution  of  the  lower  two  of  these  four 
modes  is  almost  of  negligible  importance.  The  rigid  body  mode, 
the  fundamental  mode,  and  perhaps  a  third  mode,  appear  to  suffice 
for  most  cases. 
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FREQUENCY  (CPS) 

WING  SHEAR  RESPONSE,  40.06  PERCENT  SEMIS! >AN 

Figure  B-2.  Illustrative  Frequency  Response  Function  for  Wing 
Shear  for  KC-135  Airplane  (from  reference  7) 


Firure  b-4.  Illustrative  Frequency  Response  Function  for  B-57 
Airplane  (from  reference  6) 
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WINC  KNOING  MOMENT  SPECTRUM,  2T  PERCENT  SEMI  SPAN 

Figure  B-5.  Output  Spectrum  Associated  with  Fig.  B-l  (from 
reference  7) 
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Figure  B-6.  Output  Spectrum  Associated  with  Fig.  B-2  (from 
reference  7) 
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